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ABSTRACT 


This  paper  Intends  to  be  an  Introduction  to  search  theory, 
and  It  Is  written  primarily  to  assist  staff  officers  In  their 
duties . 

For  those  who  are  in  charge  of  planning  and  deciding 
about  search.  It  Is  generally  felt  that  Is  Is  necessary  to 
know  what  Is  behind  the  standard  procedures  found  In  tactical 
publications.  However,  they  are  seldom  familiar  with  the 
advanced  mathematics  required  to  read  the  sophisticated  papers 
on  the  subject.  This  text  Is  an  attempt  to  conciliate  those 
officers  needs,  with  the  reasonable  amount  of  precision,  which 
is  required  for  planning,  analyzing  and  evaluating  search 
procedures. 
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I.  SETTING  THE  PROBLEM 


A.  SEARCH  AS  AN  OPTIMIZATION  PROBLEM 

Almost  all  relevant  search  situations  are  characterized 
by  three  basic  features: 

(a)  uncertainty  about  the  target's  position. 

(b)  limitation  and  uncertainty  about  the  capability  of 
detection. 

(c)  scarcity  of  search  resources. 

In  (a)  what  Is  meant  Is  that  target's  position  Is  described 
by  a  probability  distribution  Inside  a  "search  area"  (denoted 
"S.A. "  from  now  on) . 

In  (b)  the  point  Is  that  all  sensors  are  limited  In  their 
ranges  and  even  Inside  the  proper  range  there  still  exists  a 
positive  probability  of  missing  the  target,  due  to  factors 
as  operator's  error,  propagation  conditions,  and  so  on. 

Finally  In  (c),  the  problem  is  that  the  resource's  amount 
which  Is  available  to  perform  a  search  Is  a  finite  one.  This 
fact  Introduces  a  "budget  constraint"  In  the  problem. 

Two  basic  measures  of  effectiveness  (designated  as  "M.O.E.'s" 
from  no  on)  can  be  used  In  search  problems.  The  first  Is  the 
probability  of  finding  the  target,  In  a  given  time,  by  using 
the  available  amount  of  resources.  The  second.  Is  the  expected 
time  to  find  the  target.  As  will  be  seen,  most  of  the  pro¬ 
cedures  which  maximize  the  former,  minimize  the  latter.  For 
the  time  being,  this  paper  will  deal  only  with  the  first. 
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At  this  point,  an  example  seems  appropriate  In  order  to 
Introduce  some  Important  concepts,  and  to  Illustrate  those 
already  mentioned.  Consider  the  S.A.  In  Figure  1. 

It  Is  obvious  that  for  the  target  there  are  only  two 
possibilities:  Either  It  Is  located  In  cell  I,  or  In  cell  II 
Defining  the  events  Ej*  {Target  Is  In  cell  1}  and  E2=  {Target 
Is  In  cell  II),  the  target's  location  distribution  Is  then: 

P(E, )  -  Pj  -  0.8 

l>(E2)  ■  p„  -  0.2 

II 

Where  E  P,  *  P.  +  PTT  ■  1 
j«I  J  1  11 


Figure  1-1 

Suppose  that  the  sensor  which  Is  In  use  works  by  discrete 
glimpses,  and  define  the  events  {Target  Is  detected  In 
cell  1}  and  D2*  {Target  Is  detected  In  cell  II}.  Let 

P{0}/E})  •  Bj 

p(d2/e2)  -  b2 

These  two  probabilities  are  assumed  known  and  usually  equal, 
and  characterize  the  efficiency  of  the  sensor.  Finally,  In 
order  to  take  In  account  the  scarcity  of  searching  resources. 
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and 


P(D,)  ■  PlDj/E,)  P(E,) 

% 

And  for  cell  II,  by  Identical  reasoning 

p(d2)  ■  p(d2/e2)  p(e2) 

Let  D«  {target  Is  detected}.  Obviously,  D  Is  the  union  of 
Dj  and  D2,  and  further  more  It  Is  an  union  of  two  disjoint 
events,  for.  If  the  target  Is  detected.  It  occurs  either  In 
cell  I,  or  In  cell  II,  but  by  no  means  In  both.  Therefore, 
from  the  axioms  of  probability  theory; 

P(D1UD2)  -  P(D)  =  PfD,)  +  P(02) 
or 

P(D)  «  P(D1/E1)  P{E, )  +  P(D2/E2)  P(E2) 

In  order  to  gain  more  Insight  In  the  problem,  suppose  the 
following  numbers  are  assigned: 

Pj  «  0.8 

P„  •  0.2 

Bj  *  Bg  ■  0.6 

G  ■  2  Glimpses 

Now,  the  problem  Is  to  know  where  to  spend  the  available 
glimpses  In  order  to  have  the  larger  possible  value  for  P(D), 
the  probability  of  detection.  The  place  for  the  first  glimpse 
Is  obvious.  Cell  I  having  the  largest  probability  Is  the 
natural  candidate.  The  general  rule  Is,  evidently,  look  for 
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suppose  that  only  6  glimpses  are  permitted.  Recalling  from 
probability  theory,  the  definition  of  conditional  probability, 
It  can  be  stated  that: 


PCD-j/E-j  ) 


P(D1n£1 ) 


Where  D^rEj  stands  for  the  Intersection  of  the  events  D1  and 
Ej.  Bu^lf  this  Intersection  Is  carefully  examined.  It  Is 
not  difficult  to  see  that  It  Is  equivalent  to  the  event  D1 . 

In  order  to  prove  It,  notice  that  Is  Included  In  E-j ,  In 
the  sense  that.  If  detection  occurs  In  cell  I,  evidently  the 
target  Is  there.  Then  0^  Implies  Ejl.  On  the  other  hand,  It 
Is  possible  that  the  detection  falls  to  find  the  target  In 
Cell  I,  and  nevertheless  It  Is  there.  Then  Ej  does  not  Imply 
r1 .  Figure  2  shows  a  Venn  diagram  which  Illustrates  this  fact. 


Due  to  the  above  considerations,  (*)  becomes: 

P(D,) 

p(Di/Ei)  s  rrrfr 
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the  target  at  the  place  where  It  Is  most  likely  to  be.  The 
best  place  for  the  second  glimpse,  however.  Is  not  so  obvious. 
There  are  two  possibilities:  Cell  I  again,  or  cell  !i. 
Therefore,  by  computing  P(D)  under  these  two  possible  courses 
of  action,  It  Is  possible  to  know  which  Is  better.  Observe, 
however,  after  the  sensor  had  looked  In  cell  I,  then  Pj  and 
Pjj  no  longer  are  equal  to  P ( E ^ )  and  P ( E2 ) *  For,  If  cell  I 
was  already  scanned  and  the  target  was  not  found  there.  It  Is 
less  likely  for  this  cell  to  contain  It,  and  more  likely  for 
cell  II.  In  the  extreme  case,  when  using  a  perfect  sensor 
(Bj*B2*1),  If  the  target  Is  not  found  In  cell  I,  then  It  Is 
not  there,  and  thus  P ( E ^ )  *  0  and  P(E2)  =  1,  whatever  values 
they  had  before  the  glimpse.  Therefore,  after  the  first 
glimpse,  It  Is  necessary  to  update  P ( E ^ )  and  P(E2)  in  order 
to  compute  P(D).  Arguing  In  the  same  way  as  before.  It  can  be 
stated: 


P^/O) 


P(E,nO) 

_ t  _ 

P{D> 


Where  D  Is  the  complement  of  D.  But,  P(E^nO)  the  probability 
of  no  detection  even  though  the  target  Is  for  sure  In  cell  I, 
Is  equal  to: 

P^nO)  *  P(E1n01) 

Because  the  only  way  for  the  event  E^nB  to  occur  Is  that  no 
detection  happens  In  cell  I. 
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And  similarly 


Hence 


P(£,/0) 


_ p(e, )  n-pgyE,)] _ 

P(E,)  £l-P(D1/E-j )]  +  P(E2)  [1-P(02/E2)J 


and,  by  the  sane  reasons: 


P(E2/D) 


_ P(E2)  [i-p(d2/e2)] _ 

PU,)  C1-P(D1/E1 )]  +  P(E2)  C1-P(D2/E2)] 


The  above  results  are  a  particular  case  of  the  well  known 
theorem  of  Bayes,  from  probability  theory.  It  Is  Interesting 
to  notice  that  the  mentioned  theorem  can  be  generalized  to 
any  number  of  cells,  and  even  to  continuous  distribution. 

The  probabilities  P(E-|/0)  and  P(E2/D)  are  called  posterior 
probabilities.  In  opposition  to  Pj  and  P j j  which  are  called 
prior  probabilities.  In  section  1-2,  this  subject  will  be 
treated  In  a  more  detailed  way. 

Now,  there  Is  a  way  to  know  where  the  second  glimpse  must 
be  placed,  for,  by  using  the  numbers  of  the  problem: 

P(E./0)  -  - (P.-.8.).  LV.O-.63 -  ,  0.616 

1  (0.8)  [1-0.6]  +  (0.2)  [1-0] 

And 

p(e2/D)  -  ■+  (0.2)  Ll-aj  *  °-38i 

Therefore,  even  after  the  first  glimpse,  the  probability  of 
cell  I  Is  still  higher  and  the  second  must  be  placed  there. 

In  accordance  to  the  general  rule  stated  before. 


The  probability  of  detection,  P,  using  this  course  of 
action.  Is  now: 

P  -  (0.616)  (0.6)  >  0.3696 

Finally,  to  compute  the  probability  of  detection  on  either 
glimpse,  P(D),  by  using  the  first  option,  we  can  proceed  as 
follows: 

P(D)  ■  P  (Detection  on  I**1  glimpse)  + 

P  (no  detection  on  1th  glimpse  and  detection 
on  the  2th) 

Or 

P(D)  -  PjB,  +  PfE^O)  (1  -PjB1  ]  Bj 

0 

-  (0.8)  [(0.8)  +  (0.616)  (0.52)]  =  0.672 

Until  now,  many  Important  factors  have  been  neglected  for 
the  sake  of  simplicity.  For  example,  just  two  cells  were  used, 
no  considerations  were  made  about  the  size  of  the  search  area, 
the  sensor  was  supposed  to  work  by  separate  glimpses,  and  no 
attention  was  paid  to  the  pattern  of  the  track  that  the  search¬ 
ing  vehicle  eventually  describes  Inside  of  the  search  area. 

In  order  to  Improve  the  model,  consider  the  following  modifica¬ 
tions  In  this  example: 

Suppose  there  Is  a  target  In  the  S.A.  pictured  In  Figure 
3,  whose  position  Is  described  by  the  following  distribution: 
Probability  of  cell  j  contain  the  target  »  P j ,  jc  (I,  II,  II,  IV) 

jJj  Pj  “  PI  +  PII  +  PIII  *  PIV  “  1 
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/  \ 

*  wp 

V  ^ 

CELL  I 

CELL  II 

CELL  III 

ULL  IV  1 

Figure  1-3 


Furthermore,  suppose  that  the  sensor's  range  Is  "W"  and 
that  within  nWH  It  detects  the  target  with  probability  1. 

Let  "v"  be  the  speed  of  the  platform  where  the  sensor  is,  and 
suppose  there  Is  a  fixed  amount  of  time,  say  "T",  for  con¬ 
ducting  the  search.  If  the  area  that  can  be  swept  in  the 
time  MTM  Is  smaller  thart  "S.A.",  then  it  Is  clear  that  the 
time  must  be  shared  In  some  way  among  the  four  cells.  Of 
course,  there  are  Infinite  ways  of  doing  that,  but  the  main 
Interest  Is  the  one  that  maximize  the  probability  of  finding 
the  target.  Denoting  by  Zj  a  measure  of  effort  spent  In  cell 
"j",  and  by  MaM  the  total  amount  of  effort  available,  (In 
this  case  the  fraction  of  the  "S.A."  that  can  be  covered, 
given  "w",  Mv",  and  "T"),  It  Is  possible  to  see  that  the 
"M.O.t."  Is  a  function  of  Z  only.  Now,  the  problem  can  be 
seen  as: 
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i 


(I-I) 


MAX  P  (Zx, 


) 


S.t 


a 


Where  P(Zj,  Zjj,  Zjjj,  Zjy)  ■  P ( 2 >  Is  the  probability  of 

finding  the  target.  Zj  Is  the  DENSITY  OF  EFFORT  In  cell  j, 

l.e.,  the  ratio  between  the  effort  placed  In  cell  j,  and  the 
area  of  cell  j.  If  Aj  Is  the  area  of  cell  j  and  tj  the  time 

spent  there.  It  can  be  written  as: 


The  reason  the  density  of  effort  was  chosen  as  measure. 

Is  that  In  some  problems  there  are  cells  with  different  areas. 
Another  point  Is  that,  although  It  has  been  said  that  the 
probability  of  finding  the  target  Is  a  function  of  the  dis¬ 
tribution  of  effort  among  the  cells,  and  that  the  probability 
of  detection  Is  a  function  of  the  density  of  effort  that  has 
been  placed  there.  It  was  not  specified  any  analytical 
expression  for  It.  Coming  back  to  the  example  given,  In  order 
to  develop  an  approach  to  this  problem: 

Let  B j ( Z j )  be  the  probability  of  finding  the  target  In 
cell  j,  given  that  the  target  Is  there,  and  a  density  of  effort 
Zj  was  applied  there.  From  elementary  probability  theory 

IV 

P(Z)  ■  E  BjUj)  Pj  (I-II) 

j  *  » 


♦ 


Then,  as  Pj  Is  already  known,  the  problem  now  Is  to  find  Bj 
(Zj).  Incidentally,  Bj  (Zj)  which  stands  for  the  probability 
of  finding  the  target,  given  It  Is  there,  and  by  the  use  of 
a  density  of  effort  Zj  Is  called  a  LOCAL  DETECTION  FUNCTION, 
or  as  It  will  be  used  now,  L.D.F.  Suppose  now,  that  the 
search  Is  In  one  of  the  cells  and  assuming  that: 

(a)  If  the  target  Is  In  that  cell,  the  probability  distri¬ 
bution  of  Its  position  Is  uniform. 

(b)  If  passing  by  the  target  at  a  distance  equal  or  smaller 

than  MW"  It  will  be  detected  with  probability  1. 

7 

Under  these  circumstances,  It  seems  reasonable  to  perform  a 
systematic  search  In  the  cell.  It  can  be  done  by  using 

U 

/ 

parallel  and  non  overlapping  sweeps,  as  shown  In  Figure  I. 


f  V 

✓  * 

<  V 

Figure  1-4 


Evidently,  If  Is  the  area  of  the  cell,  the  probability 
of  detection  after  the  searching  vehicle  has  traveled  a  track 
length  "L"  Is  or  using  jp  as  the  density  effort  Z^.  (We 
will  now  skip  the  subscripts  on  Bj(Zj)): 


\WZ  If  LW<A< 
B(Z  )  j 

(l  If  LW>Aj 


(I-III) 
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This  Is  a  very  appealing  L.D.F.,  because  of  Its  simplicity. 
However,  It  Is  not  a  simple  matter  to  place  the  paths  exactly 
as  In  Figure  4. 

Even  If  It  Is  assumed  that  they  are  Indeed  parallel  one 
to  another,  errors  In  navigation  will  cause  the  paths  to  over¬ 
lap  In  some  places,  and  leave  gaps  In  others.  Thus,  the 
probability  of  finding  the  target  will  be  less  than  (I-III). 
Therefore,  unless  the  navigation  is  very  well  made,  the  use 
of  formula  (I-III)  will  overestimate  the  detection  capability. 

On  the  other  hand,  an  entire  asystematlc  procedure  can  be 
used  to  place  the  path  within  the  cell.  Instead  of  parallel 
paths,  random  ones.  It  can  be  formalized  by  adding  one  more 
assumption  to  the  other  ones: 

(c)  The  path  Is  composed  of  segments  with  random  lengths 
and  directions,  mutually  Independent. 

The  situation  Is  pictured  In  Figure  5 


The  probability  of  detection  along  segment  "AL"  Is 


f  -  (I-IV) 

Where  A  Is  the  area  of  the  cell.  Note  that  the  ratio  jp  Is 
the  effort  density  "AZ"  and  therefore  (I-IV)  can  be  written 

as 


f  -  WAZ  ( I-V) 

Let  (8 ( Z)  be  the  probability  the  target  has  been  detected 
by  the  time  the  search  vehicle  has  traveled  a  track  length 
"L",  and  therefore  placed  a  density  of  effort  "Z"  Inside  the 
cell.  Then 

B(Z  +  AZ)  -  B(Z )  +  [1  -B(Z)]W  AZ  (I-VI) 

This  Is  because,  either  the  target  Is  detected  during  the 
length  "L"  (and  B(Z)  stands  for  this  possibility),  or  the 
search  falls  to  detect  by  the  track  length  "L“  (and  succeeds 
In  the  next  Increment  AL)  and  the  probability  of'  this  event 
is  [1 -B(Z)]WAZ  (multiplication  Is  valid  here  because  of 
assumption  (c)). 

Now  rearranging  (I-VI) 

B(Z  +  AZ)  -  B(Z)  •  Cl-B(Z)]  WAZ 
or 

B.(.l+  AD  ,,B (,*.?■  «  [1  -B(Z  )]W 
tl 
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And  In  the  limit 


-  B(z)  -  B  (2)  *  [1-B(Z)]W 

/ 

Solving  this  differential  equation 

B(Z )  -  1  -  EXPC-WZ]  (I-VII) 

This  Is  called  RANDOM  SEARCH  FORMULA,  a  very  Important  one, 
first  derived  by  B.  0.  Koopman  In  [1]. 

Plotting  the  two  L.D.F's  known  already.  It  can  be  seen 
that  (I-VII)  gives  a  lower  bound  of  effectiveness  of  the  search 
Inside  the  cell.  The  graph  Is  In  Figure  1-6. 


B(2) 


The  real  L.D.F.  Is  generally  somewhere  In  the  dashed  region 
between  the  two  curves.  In  order  to  Improye  the  Insight  about 
how  precisely  either  model  describes  the  situation,  assume 
that  the  error  In  the  paths  placement  In  Figure  IV  Is  normally 
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distributed  with  mean  zero  and  standard  deviation  and  gap 
effects  mentioned  before  will  take  place,  even  though  the  paths 
remain  parallel.  It  can  be  shown  (see  [1]  and  [2]),  that,  If 
the  rectangle  Is  large  compared  to  W,  the  probability  of 
detecting  the  target  depends  only  on  the  ratio  aN/W.  Figure  7 
which  was  taken  from  (2),  shows  this  dependence. 


Ou' 

Observe  that,  for  small  ratios  y2-,  the  probability 

approaches  one  and  therefore  the  systematic  search  model  applies, 

5  N 

however,  for  the  values  of  y2-  about  0.2  the  probability 
approaches  very  rapidly  1-EXP[-1]  which  Is  the  probability 
that  would  be  found  If  WL  *  Aj  Is  substituted  In  the  formula 
( I -VI II). 

\ 

For  an  aircraft.  It  Is  normally  assumed  (see  [3])  equal 
to  20  miles  and  NW"  at  5.000  feet  of  altitude,  equal  to  40 

Ou 

miles,  which  gives  jp  ■  .5  and  the  random  search  model  applies. 


The  formula  (I-VII),  being  a  very  conservative  estimate  of 
effectiveness.  Is  far  more  useful  than  (I-III),  and  It  will 
be  used  troughout  this  paper  as  the  L.D.F.  unless  otherwise 
specified. 

Now  the  maximization  problem  (I-I)  can  be  written  In  a 
more  precise  way  because,  by  substituting  {fall)  In  (I-I)  a 
closed  expression  for  the  probability  of  detecting  the  target 
Is  obtained. 

IV 

P(Z)  »  t  P,  [1 -EXP ( -WZ  . )] 
j-I  J  J 

and  (1-1)  becomes 
IV 

MAX  l  P,  Cl-EXP ( -WZ , ) 
j-I  J  J 

s.t. 

IV 

I  A.Z.  <  a 
j-I  J  J  ~ 

To  substitute  some  numbers  In  the  example  given,  for  sake 
of  further  classification,  let:  a)  the  probability  distribution 
be: 

pi  ■  pn  '  piii  *  piv  '  T 

b)  the  searching  vehicle  speed  by: 

V  »  150  knots 

c)  the  sensor's  range  be: 

W  -  15  miles 


(I-VIII) 


(I-IX) 
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d)  the  four  cells  have  the  same  dimensions,  say  30  X  50  miles, 
and  f)  suppose  that  the  maximum  time  to  be  spent  In  the  search 
Is  no  more  than  two  hours. 


It  seems  reasonable  to  split  evenly  the  effort  among  the 
cells  and  It  Is  Indifferent  which  one  It  Is  chosen  to  start. 
The  density  of  effort  In  each  cell  Is: 

7  .  L  ,  Vt  .  100  x  i  ,1 

Zj  SJ  TO  nw1  w 

Where  t  was  used,  because  there  are  only  two  hours  to  be 
spent  In  the  search  and  y  hour  In  each  cell.  Observe  that, 
due  to  the  fact  that,  the  probabilities  and  thfe  areas  are  the 
same  for  all  the  cells  and  the  effort  Is  being  split  evenly 
among  them,  the  density  of  effort  is  the  same  In  the  whole 
search  area.  So,  It  can  be  computed  as: 

zj  JT  7TOTO  TO 

Where  the  total  time  (t*2  hours)  was  used  and  the  total  area 
(4  X  30  X  50  ■  6000  sq.  mi.).  3y  using  this  density  of  effort: 


.  1  IV  ir  l 

P(Z)  *  j  l  [1-EXP(-^)]  *  1  -  EXP[-£]  *  0.394 
*  j*I 


Now,  suppose  that  the  probability  distribution  Is  given  by 

10.4  1j  J  -  I 

0.3  1j  j  -  II 

0.2  1j  j  =  III 

0.1  1j  j  «  IV 


*  i 


And  recall  that  the  maximum  effort  Is  already  known.  The 
available  effort  Is  ■  200.  Then,  (I-IX)  becomes: 


IV 

MAX  Z  P,  [1-EXP(-WZJ 
j-I  J  J 


s.t. 

IV 

Z  A.  z.  <  200 

j-I  J  J  ” 


The  next  point  will  be  to  consider  several  ways  to  alio 
eating  the  effort  or  several  functions  Zj ,  and  see  what 
happens  to  the  probability  of  finding  the  target: 

A.  To  divide  evenly  among  the  cells  or 

zj  *  35  •  Je  ”•  n*  u,»  m 
And 

P(2)  -  I  Pj  [1-EXPO  jf)]  • 

0.4  [1 -EXP { -  ^|)]  t  0.3  [1 -E XP { -  + 

0.2  [1-EXP(-  J|)]  +  0. 1 [1 -EXP{ -  |§)]  - 
1  — EXP ( -  j)  ■  0.394 


B.  To  spend  all  effort  In  cell  I  or 


100  X  2  _  2  4 # 

W  T5  *  * :  J  I 


0  otherwise 


and 


P(Z)  -  0.4  [1  -E XP ( -  2-{s15  )]  «  0.346 

C.  To  split  evently  the  effort  between  the  cells  I  and 
II  or: 

C  ioo  x  i  _  i  ,  ft  TT, 

vt  I  ■'rare —  t?  1f  Je  {I>  111 

*  F  "S 

J  Aj  )  0  otherwise 

and 

P(Z)  ■  0.4  [1  -EXP( -  ||)]  +  0.3  [1  -EXP ( -  j|)]. 

-  0.4431 


D.  To  divide  the  effort  proportionally  to  the  probability 


of  each  cell  or 


0.0533  If  j  *  I 

0.0399  If  j  *  II 

0.0266  If  j  *  III 

0.0133  If  j  *  IV 


Where  1.66  Is  the  speed  In  miles  per  minute  and  48,  36,  24, 

12  the  minutes  spent  In  cells  I  to  IV  respectively.  The  pro¬ 
bability  Is: 


P(Z)  =  (0.4)  (1 -EXP ( -0 .8)3  +  (0.3)  [l-EXP(-0.6)]  + 

(0.2)  [1 -EXP (-0.4)3  +  (0.1)  [1 -EXP ( —0 . 2 ) ]  = 
0.4402 
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E. 


To  divide  evenly  the  effort  amongthe  cells  I,  II  and 


III  or 


7  .  vt 

zj  *7 


•4'°'T5M66  s  °-0442  1f  JEU.  ”•  III} 
0  otherwise 


and  by  similar  calculations 


P(Z)  =  0.4356 


F.  To  divide  the  effort  proportionally  to  probability, 
among  cells  I,  II  and  III;  do  not  place  any  effort  in  cell  IV 
or 

0.058  If  j*I 

0.0442  If.  II  • 

0.030  If  J*III 
^  0  If  j  -  IV. 
and  the  probability  Is 

P(z)  *  0.450 

The  table  In  Figure  8  summarizes  the  results  reached  till 
now,  by  the  use  of  six  courses  of  action. 


if 


Vt 

*7 


53  X  1.66 
“TSffi - 

40  X  1.66 
U00 

27  X  1.66 

~rm — 


ALLOCATION 


PROBABILITY  OF 
FINDING  THE  TARGET 


A  P(A)  •  0.3940 

B  P(B)  -  0.3460 

C  P(C)  -  0.4431 

D  P(D)  -  0.4402 

E  P(E)  -  0.4356 

F  P(F)  -  0.4500 

Figure  1-8 

It  can  be  seen  that  for  different  allocations,  the  pro¬ 
bability  varies  In  a  wide  range,  (the  range  could  be  wider 
If  some  less  obvious  allocation  was  chosen,  say,  to  spend  all 
time  In  cell  IV),  and  that  common  sense  by  Itself  Is  unable 
to  provide  any  guidance,  even  In  the  simple  example  we  are 
using.  For,  the  best  result  achieved  was  by  placing  no  effort 
In  cell  IV,  which  does  not  seem  very  wise  at  a  first  sight. 
Furthermore,  there  Is  not  any  procedure  to  know  If  we  can  do 
any  better,  than  the  one  done  with  allocation  F,  because 
there  are  Infinitely  many  possible  allocations  and  we  have 
only  checked  six. 

In  the  next  section,  a  specific  procedure  will  be  considered, 
which  works  very  well  In  most  cases. 

Another  point  seen  In  the  example,  Is  that  the  same 
amount  of  effort  (200)  has  been  used  In  all  six  allocations. 


It  can  be  said  that  all  of  them  "cost"  the  same,  but  "F" 
gives  a  better  NreturnN  than  the  others.  In  fact.  It  can  be 
defl ned 


J 

E  A,  Z,  -  C  (Z  ) 
J-I  J  J 


As  the  cost  of  allocation  Zj ,  and  C  (Z)  will  be  used  to 
denote  the  function  and  "CN  to  denote  the  values  It  takes. 
Then,  generalizing  (1-1),  the  search  problem  can  be  set  as: 

MAX  P(Z) 


s.t. 

C  (2)  <  C 


(I-X) 


Or  putting  It  In  words,  on  allocation  Is  desired,  that  Is 
a  solution  for  the  optimization  under  constraint  problem 
(I-X). 


B.  POSTERIOR  TARGET'S  LOCATION  DISTRIBUTION 

Returning  to  the  example  given.  It  can  be  recalled  that 
Pj  was  defined  as  the  probability  that  the  target  Is  located 
In  cell  j.  Howevar,  these  probabilities  are  computed  before 
any  search  effort  Is  placed  In  the  S.A.  Once,  some  effort  Is 
placed,  those  values  no  longer  hold,  for  after  each  Increment 
of  effort  placed  In  a  cell  without  success,  the  probability 
that  the  target  Is  In  the  searched  cell  decreases,  and,  by 
consequence,  the  other  cells  probabilities  will  be  Increased. 
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Then,  Pj  can  be  said  as  the  probability  associated  with  cell 
j  after  some  effort  was  placed  In  the  S.A.  "P  Is  called  the 
prior  probability  and  "PjM  the  posterior  probability,  and  the 
corresponding  distributions  prior  and  posterior  target  loca¬ 
tion  distribution,  respectively. 

From  probability  theory,  Bayes'  theorem  Is  recalled,  (for 
demonstration  see  [4]),  and  by  Its  use  posterior  probability 
Is  computed  as  follows:  Suppose  a  search  has  been  made  using 
allocation  and  that  we  did  not  find  the  target.  Then: 


Pk  [l-B(Zk)J 

pk  '  3 - 

i  p,  n-BU,)] 
j-i  3  3 


(I-XI) 


Where  the  numerator  stands  for  the  probability  that  the  target 
Is  located  In  cell  "K",  and  Is  not  detected,  and  the  denomina¬ 
tor  the  probability  of  no  detection  In  any  cell. 

For  example,  suppose  a  search  was  made  according  to 
allocation  "B"  and  failed  to  find  the  target  after  the  two 
hours  established.  The  posterior  distribution  Is  computed: 

Pk  [EXP(-ZkW)] 

’’k  '  TV - - - 

r  Pi  [exp(-z,w)] 
j-i  J  3 

Then: 

•  *  .  (0,4)  cexp<-  nrffii 

1  (0.4)  C EXP ( -  +  0.6 


0.08 


0.46 


;  .  _ Q.3  _ 

11  (0.4/  [EXP(-^^-)]  +  oTs 


III 


0.2 


(0.4)  [EXP( 


Trrrrr 

t r'J 


+  0.6 


0.310 


Piv  "  - HtTT -  *  0.153 

IV  (0.4)  [EXP ( -  ^Y‘)  +  0.6 

The  table  below  displays  the  effect  on  the  target  position 
distribution  In  two  hours  of  search  In  cell  I 


BEFORE 

AFTER 

j  ■  o-i 

Pj  «  0.080 

II  *  0,3 

Pn  -  0.460 

III  "  °*2 

Pin  -  0.310 

IV  "  °*1 

Pjy  -  0.150 

Almost  all  the  probability  mass  which  was  in  cell  I  in 
the  very  beginning,  was  transfered  to  the  other  cells.  There¬ 
fore,  If  there  Is  more  time  left  to  continue  the  search,  It  Is 
reasonable  to  transfer  the  effort  from  cell  I,  because  It  Is 
almost  Impossible  the  target  Is  located  there.  Indeed,  a  very 
Intuitive  way  of  allocated  effort  Is  to  place  It  where  the 
probability  of  success  Is  higher.  If  the  search  starts  In  the 
cell  with  highest  prior  probability,  and  the  criterion  of 
placing  the  next  Increment  of  effort  In  the  cell  with  highest 
posterior  probability  Is  used,  an  allocation  Is  defined.  As 
will  be  seen  later,  this  allocation  Is  In  most  cases,  (Including 
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Then, 


In  this  case 


A1C1  _  30X50X0.01918 


0.28770  hours 


The  posterior  distribution,  arter  e ^ : 
0.4  EXP( -1 5e-, ) 

fi  _  "  I 


T  “  7 


uS  "  I 


PIV  • 


nr  *  9 


Now,  In  accordance  to  the  procedure  chosen,  the  effort 
must  be  split  evenly  between  cells  I  and  II  until  another  cell 
(In  this  example  cell  III)  shows  a  higher  posterior  probability. 
Then,  for  computing  e2: 


pi  ■  pm 


EXP  [-15e2] 


2(575)  EXP  [-15c2]  +  f  +  5 


Where  P,  stands  for  the  posterior  probability  of  cell  "j 

«** 

after  e2  has  been  placed.  Equating  Pj  or  P  j  j  to  Pjjj, 


(^-gO  EXP  [ - 1 5 Eg 3  *  g- 


rCc  e  ,  c 3  i  1 


because,  as  always,  the  denominators  are  Identical.  Then, 


1n  (f) 

15 


0.027 


And  the  time  Increment 

’2 7  e2  .  2X30X50X0.027 


0, 


(A^  +  A9)  e. 


100 


«  0.81  hours 


To  be  spent  half  in  cell  I  and  half  In  cell  II.  The  new 
posterior  distribution  after  the  density  e2  Is: 


PT  «  P 


II  rIII 


2  (jt|)  EXP  [-15e2]  ♦§  +  \ 


(ff7§)  EXP  [-15e2] 


2 

7 


:  1 

U  ’  2(§t|-)  EXP  [-15e2]  +  5  ’  7 


In  the  same  way  we  compute  e3,  which  Is  the  effort  density  to 
be  placed  In  cells  I,  II  and  III  before  moving  to  cell  IV. 


pi  ■  pn  ■  piii 


(f)  EXP  [-15e3] 

3(f  EXP  [-15e3]  +  } 


(f)  EXP  [-15e3]  -  } 


and 


In  (j-) 


15 


0.04621 
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Observe  that,  In  the  example  given,  there  Is  a  limit  of 
two  hours  to  perform  the  search.  Then,  all  0j  and  0^  can  be 
used,  but  only  a  fraction  of  03»  that  will  be  called  0^.  It 
Is  obvious  that: 

-  2  -  (01  ♦  02)  »  2  -  (0.28770  -  0.81000) 

■  0.90230  hours 

So,  the  time  to  be  spent  In  each  cell  will  be 

0  0* 

1 1  “  0 1  +  j-  B  0 .99347  hours 

02  ei 

lz  a  T'  +  T  9  °-70576 
6 1 

t3  -  y-  -  0.30076 


Where  "tj"  stand  for  the  total  time  In  cell  Mj".  Observe 
also  that  t ^ *Z j  in  the  example  given,  due  to  the  numbers 
chosen.  Hence  the  probability  of  finding  the  target  can  be 
computed  as 

P(2)  »  (0.4)  [1  -EXP ( -t] ) ]  +  (0.3)  [l-EXP(-t2)]  + 

(0.2)  [l-EXP(-t3)3  +  (0.1)  [1 -EXPf-t^)] 

-  0.4553 
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Which  Is  better  than  any  other  of  six  allocation  tried  In  the 
last  section.  In  fact,  as  said  before.  It  Is  the  best  that 
can  be  done  under  these  circumstances.  Finally,  It  Is 
Interesting  to  observe  that  using  the  procedure  above.  If  for 
any  reason  there  Is  a  need  to  stop  before  all  the  Initial 
effort  available  Is  gone,  the  M.O.E.  has  still  been  maximized, 
under  the  new  and  smaller  constraint.  On  the  other  hand.  If 
at  the  end  of  the  operation  more  effort  Is  available,  then 
the  procedure  chosen  need  no  change;  the  search  continues  by- 
placing  the  effort  where  the  posterior  probability  Is  highest. 
Therefore,  the  procedure  has  a  sort  of  Independence  In  relation 
the  constraint.  This  point  Is  a  very  Important  feature  in 
real  world's  operations  (a  change  in  the  weather,  as  example, 
can  determine  the  end  of  search)  and  characterizes  what  will 
be  defined  later  as  an  UNIFORMLY  OPTIMAL  SEARCH  PLAN. 

C.  SWEEP  WIDTH 

In  the  preceding  section,  a  quantity  "W"  equal  to  15  miles, 
was  used  as  a  distance  within  which  the  sensor  will  detect  a 
target  with  probability  1.  In  the  real  world  however,  things 
are  not  so  simple.  From  a  rigorous  point  of  view,  the  detec¬ 
tion  ability  of  a  sensor  can  only  be  described  by  a  random 
variable  which  associates  to  each  distance  a  probability  of 
detection,  for  a  given  type  of  target,  a  given  condition  of 
propagation  (light,  sound,  Infrared,  radar  waves,  etc.)  a  given 
level  of  personnel's  aptitude  and  a  multitude  of  other,  perhaps 
less  Important,  factors.  It  should  be  obvious  that  It  Is  not 
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and 


(Aj  +  ^2  +  A3)  c3 
V 

3XJL0.X50X_0,0462J,  .  hours 

100 


To  be  spent  one  third  In  each  one  of  the  cells  I.  II,  and  III 
The  posterior  distribution  after  e3  Is: 


pn  *  pm 


(7)  EXP  C-15e3] 

4  (f)  EXP  [-15  e3] 


1 

*  ‘ 


Figure  9  displays  the  variation  of  the  cell's  probabilities, 


Before  the  After  e. 

Search 


Figure  1-9 
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easy  to  obtain  the  distribution  of  such  a  random  variable. 

In  fact,  In  some  attempts  to  derive  the  probability  of  detec¬ 
tion  for  search  radar,  considering  a  given  distance,  the 
variance  of  the  sample  was  so  big,  that  the  corresponding 
confidence  Internal  for  the  estimate  was  wider  than  would 
reasonably  be  acceptable.  The  way  normally  used  to  overcome 
such  difficulties,  Is  by  means  of  a  lateral  range  curve,  from 
which  a  SWEEP  WIDTH  Is  derived. 

Consider  Figure  10  where  the  vertical  lines  represent  the 
target's  paths.  The  distance  between  the  sensor  and  the 
Intersection  of  a  vertical  line  with  the  horizontal  axis,  is 
the  lateral  range  of  the  corresponding  path. 


Pi 

N 


Figure  I - 1 0 
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If  several  targets  are  passed  through  each  path  and  the 
percentage  of  those  detected  computed,  no  matter  In  what  point 
of  the  path  It  happens,  this  percentage  can  be  taken  as  the 


probability  of  detection  of  a  target  for  which  the  closest 
point  of  approach  is  that  lateral  range.  By  plotting  those 
probabilities  against  distance,  the  lateral  range  curve  for 
the  sensor  In  question  Is  obtained.  An  example  of  lateral 
range  curve  for  a  radar.  Figure  1-11  shows: 


P(*) 


X 


Figure  1-11 

Observe  that  the  lateral  range  curve  does  not  represent 
a  probability  distribution,  observe  also  that  the  probability 
drops  near  the  origin.  This  Is  due  to  the  sea  return. 
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The  SWEEP  WIDTH,  which  will  be  denoted  by  W  Is  the  area 

below  the  lateral  range  curve.  The  "W"  used  In  sections  I  and 

II  corresponds  to  a  lateral  range  curve  as  the  one  In  Figure 

w 

12,  where  R  »  £  is  called  the  range  of  a  definite  range  law 
of  detection,  a  simplification  very  much  used  In  real 
situations . 
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Figure  1-12 

Although  In  chapter  II  this  subject  will  be  discussed  again. 
It  is  Interesting  to  remark  now,  that  In  some  situations, 
(chiefly  In  under-water  search).  It  is  necessary  to  treat  the 
sweep  width  as  a  random  variable. 
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0.  MOVING  TARGETS 

Until  now,  all  the  derivations  and  concepts  presented  In 
this  paper,  were  developed  under  the  assumption  that  the 
target  remains  stationary  during  the  search.  If  It  Is  wanted 
to  consider  moving  targets,  the  following  remarks  are  Important: 

(a)  The  problem  of  effort  allocation  under  these  circum¬ 
stances,  Is  a  very  difficult  one,  and  for  most  of  the  cases, 
has  not  been  solved  yet.  The  reason  for  this  Is  that  the 
prior  probability  distribution  for  the  target's  location  Is 
modified  continuously  In  time,  In  an  unknown  way,  by  the 
target's  movements.  A  case  where  the  problem  Is  tractable 
Is  where  the  target  motion  Is  deterministic,  and  with  known 
parameters.  In  this  case,  the  target  "drags"  the  probability 
distribution  In  a  known  way,  and  as  will  be  seen  later.  It  Is 
possible  to  solve  the  effort  allocation  problem. 

(b)  If,  however,  the  search  Is  to  be  made  at  random,  then, 
the  preceding  discussion  for  the  random  search  applies,  what¬ 
ever  prior  distribution  Is  assigned  to  the  target's  location, 
and  for  any  kind  of  target's  motion.  The  rational  Is  that. 

If  the  sensor's  movements  are  at  random  then  any  fraction  of 
the  search  area  Is  equally  likely  to  be  scanned  at  a  given 
time,  and  evidently  It  applies  to  the  place  where  that  target 
Is  at  the  considered  time. 
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E.  THE  CONTINUOUS  CASE 

Suppose  that  Instead  of  cells  the  S.A.  Is  thought  as  a 
continuous  set  of  points.  Then,  Instead  of  cells'  probabili¬ 
ties  we  will  deal  with  points'  probabilities,  or  more  exactly, 
with  a  continuous  bldlmenslonal  probability  distribution. 

Observe  that  In  this  case,  most  of  our  definitions  should 
be  reformulated.  An  allocation,  for  example,  will  be  a  function 
which  assigns  to  each  point  of  the  S.A.  an  effort  density.  In 
a  similar  way  L.D.F.,  and  so  on,  must  be  redefined.  This  will 
be  done  later,  but  for  the  time  being  two  points  must  be 
emphasized: 

First,  for  theoretical  purposes,  the  continuous  approach 
Is  better,  not  only  because  It  Is  closer  to  the  real  situation 
(after  all  a  S.A.  Is  continuous  set  of  points),  but  also 
because  of  the  mathematical  tractablllty  of  continuous  function, 
which  permits  us  to  derive  general  results  and  theorems. 

Second,  for  practical  purposes,  the  continuum  of  the  S.A. 

Is  approximated  by  a  large  number  of  small  cells  and  then  the 
theory  developed  through  the  continuous  model  approximately 
applies. 

Hence,  all  the  Intuitive  discussion  of  I-A  to  I-D  are 
entirely  justified  and  the  more  sophisticated  discussion  In 
following  pages  will  have  the  sole  purpose  of  Improving 
foundations . 
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II.  BASIC  THEORY 


Under  this  heading,  what  Is  aimed  Is  not  only  to 
formalize  from  a  mathematical  point  of  view  the  Intuitive 
notions  Introduced  In  I,  but  also  to  provide  u  general 
overview  of  the  main  results  of  search  theory. 

A.  SEARCH  SPACE 

As  was  mentioned  before,  the  target's  location  Is 
described  by  a  probability  distribution.  More  precisely, 

It  Is  said  that  this  location  Is  a  random  variable. 

The  range  of  this  random  variable,  (or  better,  the 
essential  range,  see  [4]  ),  which  Is  some  subset  of  the 
plane.  Is  called  the  search  space. 

Observe  that  the  search  space  can  be  Infinite,  (as  will 
be  shown  for  the  case  of  the  bldlmensl onal  normal  distribu¬ 
tion,  for  example),  and  this  is  not  a  convenient  feature,  at 
least  In  practical  situations.  Because  of  this,  search 
area  Is  defined  as: 

B.  SEARCH  AREA 

Is  the  subset  of  the  search  space  where,  search  operations 
will  be  conducted. 

As  was  shown  In  definition  II-A,  the  search  space,  being 
any  subset  of  the  plane,  can  In  fact  be  very  pathological. 
Consider,  as  an  example,  the  search  space  represented  In 
Figure  1 1 - 1 . 
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Figure  1 1 - 1 


The  concentric  circles  stand  for  Isoprobablllty  curves. 

It  Is  a  case  of  a  Infinite  search  space.  If  the  dashed 
area  encloses  a  satisfactory  amount  of  probability  for 
the  purposes  sought,  the  operations  can  be  confined  to  that 
same  area.  The  dashed  portion  of  Figure  1  Is  then  an 
example  of  search  area. 

Finally,  In  accordance  with  the  model  used,  the  search 
area  can  be  dealt  with  as  a  continuum,  or  It  can  be  divided 
in  an  arbitrary  number  of  cells,  assigning  to  each  one  the 
probability  mass  It  encloses.  Furthermore  If  cell  "k" 
encloses  probability  It  Is  assumed  that  this  amount  of 
probability  Is  uniformly  distributed  In  the  cell. 

C.  LOCAL  DETECTION  FUNCTION 

Let  "x"  be  a  point  In  the  search  area  and  "Z."  the  effort 
density  placed  there.  The  function  which  to  each  pair  (x,Z) 
assigns  the  probability  of  detecting  a  target  located  In  x  Is 
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called  a  local  detection  function.  If  points  are  replaced 
by  cells  and  for  sake  of  notatlonal  consistency  "x"  by  "j" 

In  the  preceding  lines*  the  definition  given  holds  for  the 
descrete  case. 

Note  that  "x"  Is  a  bldlmenslonal  vector ,  x  *(x^ ,  x2)  and 
"2"  a  non-negative  real  number.  The  L.D.F.  will  be  referred 
to  as  B(x*Z).  This  notation  takes  In  account  the  possibility 
(very  remote  In  real  world  problems)  of  different  forms  of 
L.D.F.  In  different  cells. 

D.  L.D.F.S  AND  DETECTION  RATES 

In  I  two  L.D.F.s  were  derived  one  for  systematic  and 
the  other  for  random  search.  Here*  a  more  precise  approach 
will  be  used  to  obtain  a  general  formulation  from  which  the 
above  mentioned  L.D.F.s  are  particular  cases. 

Consider  a  detection  system  which  operation  consists  of 
a  series  of  Instantaneous  observations*  for  example,  measuring 
a  distance  to  a  submarine  by  a  step  scanning  sonar.  It  Is 
called  a  discrete  system.  The  mentioned  observations  will 
be  called  glimpse  and  let  "g"  denote  the  probability  of 
detection  on  a  glimpse.  If  "gH  Is  a  constant,  the  glimpses 
are  Independent  and  If  "N"  denotes  the  random  variable  "number 
of  glimpses  until  detection  occurs".  Then: 

P  (N*n)  *  g  (l-g)n-1  ( I I -I ) 

P  (N<n)  -  1  -  (1  -g )n  (II-II) 
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Because  HNN  Is  evidently  a  geometric  distributed  random 
variable  and  (II-I)  and  (II-II)  are  respectively  Its  probability 
mass  function  and  cumulative  distribution  function.  Suppose 
now  that  the  system  used  Is  a  continuous  one,  1,e,  It  does 
not  operate  only  at  discrete  points  In  the  time,  but  Instead 
It  may  reach  detection  at  any  moment.  Rather  than  define 
and  Instantaneous  probability  as  done  before,  let  y  be  a 
"detection  rate"  defined  in  such  a  way  that  yAt  Is  the 
probability  of  detecting  the  target  during  an  Interval  of 
length  At,  given  It  has  not  been  detected  before.  In  this 
case  let  "T"  be  a  random  variable  which  stands  for  the  time 
at  which  detection  occurs.  In  analogy  with  (II-II),  It  Is 
Intended  to  find: 

P  (t)  -  P  (T<t)  (II-III) 

consider 

Q(t)  -  l-P(t)  -  P(T>t)  (II-IV) 

the  "survival  function"  of  "T".  It  Is  clear  that  the  relation 

Q(t+At)  *  Q(t)  (1-yAt)  '  (II-V) 

holds,  by  arguing  as  follows:  Q(t+At)  is  the  probability 
of  the  event  (detection  does  not  occur  by  the  time  t+  At), 
which  Is  logically  equivalent  to  the  Intersection  of  the 
events  A=(detect1on  does  not  occur  by  time  t)  and  B=(detect1on 
does  not  occur  during  the  Interval  (t+  At  )). 

But  by  definition: 

P(A)  =  Q(t ) 

and  recalling  that  the  Interval  (t  +  At)  has  length  "At",  and 
that  "yAt"  stands  for  the  conditional  probability  that  detection 
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has  occured  during  the  Interval*  given  It  has  not  happened 
before: 

P(B/A)  «  1-yAt 
from  probability  theory: 

P(AnB)  »  P ( A )  P ( B/A ) ,  hence  (II-V) 

A  simple  manipulation  of  (II-V)  yields: 

Q(ttAt)  -  Q ( t )  .  .  Y„(t) 

or  In  the  limit  when  At  -  0 
Q'(t)  *  -  rQ(t) 

whose  solution  Is: 

Q(t)  ■  Q(0)  EXP  C-yt] 

Since 

Q  ( 0 )  -  1 
Obviously 

P(t)  -  1  -  EXP  C-yt]  (II -VI ) 

Observe,  however,  that  both  (II-II)  and  ( 1 1 - VI )  have  been 

derived  under  the  assumption  that  "g"  and  Y  remain  costant 
all  time.  If  Its  variation  with  time  Is  to  be  taken  In 
account,  the  expressions  can  be  modified  by  writing: 

P(N<n)  •  1  { 1 -g, )  (II-VII) 

for  the  discrete  case,  or 

P(t)  -  1  -  EXP  [-  f  y ( t )  dt]  (II-VIII) 

for  the  continuous,  where  "g^H  stands  for  the  detection 
probability  on  the  1^  glimpse,  and  y ( t )  Is  a  time  dependent 
detection  rate.  Although  there  are  many  reasons  by  which  g 
and  or  "y"  do  not  remain  constant  In  time  (a  human  factor 
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reason  Is  the  fact  that  the  device's  operator  "learns"  while 
he  Is  watching  the  display),  only  the  distance  variation  will 
be  considered  as  a  source  for  the  time  dependence  of  "g"  or 
"y".  Thus,  If  by  denoting  r^  as  the  distance  to  the  target 
at  the  time  of  1th  glimpse,  and  "r  (t)"  the  distance  at  time 

N 

91  ■  g(r1 ) 
and 

y( t )  -  Y[r(t)] 

By  substituting  In  ( 1 1- VI I )  and  (II-VIII) 

P(N<r)  -  1  -  If  (l-g(r,))  (II-IX) 

“  1>1  1 

and  1  1 

P(t)  -  1  -  EXP  [-  /*  y  Cr(t)]  dt]  (II-X) 

o 

notice  that  the  probability  of  detection  during  the  Interval 
(tQ ,  tQ  +dt)  can  be  computed  If  It  Is  a  product  of  probability 
of  no  detection  before  tQ  (which  Is  given  by  1  minus  formula 
(II-X),  times  probability  of  a  detection  during  "dt"  (which 
Is  Y  Cr ( tQ ) ]  dt),  or: 

Y  [r(tj]  dt  EXP  [-  /  0  >[r(t)]  dt] 

0 

but,  on  the  other  hand.  It  Is  also  equal  to  the  difference 

P(tQ  ♦  dt)  -  P(tQ)  :  P'  (tQ)  dt 

thus  f 

lo 

P'  (O  dt  -  y  [r(t  )]  dt  EXP  [-  /  y  Cr(t)]  dt] 

°  u  o 

and 

P'(tJ  ■  y  [r(t  )]  EXP  [-  /°  y  [r(t)]  dt]  (II-XI) 

0 
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P'(t0)  Is  the  value  at  tg  of  the  density  function  of  random 
variable  "time  until  detection". 

In  Figure  II- 2  ,  which  Is  taken  part  from  [1],  part  from 
[5]  five  examples  of  Y[r*(t)]  are  shown: 

A)  case  I  applies  when  (r)  reaches  a  finite  maximum  at 
distance  zero: 

B)  case  II  Is  when  the  maximum  Is  Infinite, 

C)  case  III  shows  the  effect  of  sea  return  shifting  the 
maximum  away  from  the  origin. 

D)  case  IV  Is  when  the  probability  Is  one  when  r<rQ. 

E)  case  V  Is  when  the  so  called  definite  law  of  detection 
occurs.  In  this  case,  detection  Is  sure  to  occur  as  soon  as 
the  target  reaches  the  dashed  region  and  Is  Impossible  outside 
of  it. 


Figure  1 1 -2 
45  . 
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Another  Important  form  of  detection  rate  used  for  the  case  of 
visual  detection  Is  "the  Inverse  cube  law  of  sighting", 
whose  expression  Is: 


Y  (r) 


Kh 

(hz  ♦  r2)3'2 


(II-XII) 


Where  "h"  Is  the  observer's  height,  "k"  a  constant  which 
depends  on  fixed  factors  as  contrast  of  wake  against  ocean, 
observer's  ability,  metereologlcal  conditions,  area  of  the 

target's  wakes  etc . .  and  "r"  Is  the  distance.  The  above 

formula  Is  derived  under  the  assumptions  that  the  observer 
detects  the  target  by  seeing  Its  wake,  and  "  Y(r)"  is  propor¬ 
tional  to  the  solid  angle  subtended  at  the  point  of  observation 
by  the  wake. 

When  "r"  Is  much  large  that  "h"  (II-XII)  becomes: 


T  (r)  »  (II-XIJI) 

rJ 

The  details  concerning  the  derivation  of  (11-12)  are  In  [1], 

E.  LATERAL  RANGE  CURVES  AND  SWEEP  WIDTH. 

Although  some  nice  expressions  for  y(r)were  derived  In 
last  section,  nothing  was  said  yet  about  how  to  deal  with  r, 
the  distance  between  the  target  and  the  sensor. 

It  Is  apparent,  however,  that  It  Is  very  difficult  to  keep 
track  of  the  distance  variation  when  performing  a  search.  For, 
If  there  Is  relative  motion  between  target  and  sensor,  at 
each  Instant,  r,  will  assume  a  different  value,  and  therefore 
P(t)  will  depend  In  the  law  of  variation  of  r,  as  function  of 
time.  This  Is  extremely  Inconvenient,  because  even  for  very 
simple  types  of  motion,  this  problem  can  be  an  extraordinarily 
difficult  one. 
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In  order  to  avoid  the  mentioned  difficulty  some  theory 
must  be  developed.  Suppose  there  Is  a  need  to  compute  the 
probability  of  detection  along  a  given  track.  A  rectangular 
(Xj,  x g )  coordinate  system,  will  be  adopted,  with  the  sensor 
at  the  origin  (0,0).  Notice  that  a  relative  motion  approach 
Is  being  used.  The  scalar  distance  from  the  target  to  the 
sensor.  Is  a  function  of  the  time  and  Is  denoted  ad  r  (t). 
if  the  track  that  Interest  Is  the  one  In  Figure  3,  It  Is 
easy  to  see  that  r(t)  *  /XZ  (t)  +  „Z  (t)  . 


At  this  point  It  Is  convenient  to  define  a  quantity 
"F",  which  Is  a  characteristic  of  each  possible  path,  and 
which  Is  known  as  the  sighting  potential  along  the  considered 
track.  In  order  to  make  clear,  the  dependence  of  the  sighting 
potential  on  the  track  It  will  be  denoted  by  F(.),  where  the 
dot  Is  to  be  replaced  by  some  Indicator  of  which  track  is 
being  considered  at  the  moment.  Example:  Suppose  the  paths 
are  numbered.  Then  F(1)  refers  to  path  1,  and  so  on.  By 
dlflnl tlon: 

F  (*)  *  y  1  (t>  +  (t)]  dt 

For  the  case  of  the  continuous  looking,  and 

F  (•)  *  in  [1-g  (  /x2  (t)  +  xz  (t)>] 

for  discrete  model.  By  the  use  of  the  F(.)  concept,  equations 
(II-1X)  and  ( 1 1 - X)  may  be  united  Into 

P(-)  -  1  -  EXP  [-F( • )] 

Which  gives  the  probability  of  detection  along  the  considered 
track,  no  matter  If  a  continuous  looking  Is  being  used,  or 
separate  glimpses.  The  verification  Is  a  matter  of  trivial 
algebra.  It  Is  Important  to  state  the  property  of  additivity 
which  the  sighting  potentials  have.  Suppose  the  probability  of 
detection  along  at  least  one  of  two  Independents  tracks  (say 
track  1  and  track  2)  Is  to  be  computed,  then: 

P(l)  =  1  -  EXP  C-F ( 1 )] 

P (2)  «  1  -  EXP  C-F(2)] 
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The  only  way  to  fall  In  detecting  is  to  fall  In  both  paths  or 
1  -  P(1  or  2)  =  [1-P(1 )]  [1 -P(2 ) ] 
but 

Cl-P(l)3  Cl-P(2)]  -  (EXP[-F(1 )])  (EXP  [-F(2)3)  - 

EXP  [ - ( F ( 1  )  +  F(2))3 

then  P(1  or  2)  *  1  -  EXP  C-(F(1)  +  F(2))3 

what  means  that  If  F(l)  and  F(2)  are  known,  just  by  adding 
them,  the  sighting  potential  along  the  total  path  Is  obtained. 

The  additivity  applies  to  any  number  of  paths  and  can  be 
useful  to  deal  with  complicated  tracks,  by  decomposing  them  In 
a  series  of  simpler  ones,  whose  potential  can  be  easily  computed 
and  just  by  adding  the  sighting  potentials. 

Consider  now  the  point  labelled  "CPA"  In  Figure  3.  The 
Initials  stand  for  closest  point  of  approach,  and  the  distance 
x  from  the  CPA  to  the  sensor.  Is  the  target'  lateral  range. 
Observe  that  y(r)  Is  a  function  of  the  distance  and  therefore, 
for  all  practical  purposes.  It  Is  zero,  beyond  a  certain 
distance.  Thus  It  only  makes  sense  to  consider  the  variation  of 
y(r)  for  values  of  r  between  x  and  rQ  such  that  y(r)*0  If  r>rQ. 
Suppose  that  we  replace  the  actual  path  In  figure  3,  by  Its 
tangent  at  the  CPA.  The  straight  line  being  as  long  as  we 
want,  all  distances  between  x  and  rQ  can  be  considered.  On 
the  other  hand,  the  linearity  of  the  path  makes  the  problem 
much  simpler.  For  example  to  compute  F(.),  instead  of  an 
Integral  along  a  path,  a  common  Rleman  Integral  Is  enough 
Finally,  observe  that  when  the  relative  velocity  between  sensor 
and  target  Is  constant,  and  the  most  common  example  of  this 
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case  Is  a  sensor  moving  at  constant  velocity  and  a  stationary 
target,  then  the  relative  track  Is  actually  a  straight  line. 

Consider  now  Figure  4,  where  a  straight  line  path  Is 
represented. 


Figure  II-4 

as  It  Is  already  known: 


'r(t)  ■  'x*  (t)  ♦  x\  (t) 

but  now  x1  (t)  Is  a  constant,  "x"  the  lateral  range  and  as  seen 
In  figure  4^,  Is  the  distance  of  the  closest  approach  between 
the  target  and  the  sensor.  Then, 


furthermore  as  the  direction  of  the  relative  velocity  vector 
Is  parallel  to  the  "x2”  axis  (the  axis  were  chosen  In  this 


r(‘»  *  '7777 


F(-)  •  /  r  (Q  t  v'?t"?)  dt 
F(-)  •  -  I  In  [1-g  <'7777>] 

For  the  continuous  looking  case  and  separate  glimpses  case 
respectively. 

Observe  that  now,  depending  on  the  form  of  the  detection 
rate,  both  the  Integral  and  the  sum  above  looks  like  math¬ 
ematically  tractable.  In  fact  for  many  cases  they  have  been 
evaluated  In  simple  mathematical  form.  The  Interested  reader 
can  see  [1]  for  a  complete  discussion.  For  the  sake  of 
further  simplification,  let  the  case  where  the  path  Is  an 
Infinitely  long  straight  line  be  examined.  In  this  case  t'-*--* 
and  t"-*-»  and  the  sum,  which  is  used  in  the  definition  of  the 
separate  glimpses  cases,  becomes  a  series.  In  both  cases 
however,  F(.)  Is  now  a  function  of  the  lateral  range  only. 

Thus  It  Is  more  appropriate  to  write  F(x)  Instead  of  F( . ) . 

The  probability  of  detection  along  such  a  path,  which  by  the 
same  reasons  Is  also  a  function  of  "x"  only,  can  be  computed 
P(x)  »  1  -  EXP  C-F(x)] 

which  defines  a  function  called  a  lateral  range  curve.  There¬ 
fore,  even  though  we  cannot  truly  handle  general  paths,  we  can 
approximate  them  by  use  of  P(x). 

Hence  forward,  In  this  paper  only  the  continuous  looking 
case  will  be  dealt  with.  The  Interested  reader  can  consult 
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[1]  or  [5]  for  results  and  discussion  concerned  to  the  separate 

glimpses  case.  As  It  has  been  shown  before 

t" 

F(*)  -  /t  y  (/XZ  +  yZtZ)  dt 
when  t and  t"** •  : 

f(x)  -  /"  y  (S7T7?)  dt 

or,  considering  that 
vt  ■  x2 
It  Is  obtained 


F(x)  -  7  /  Y  1^7,!')  d  x 


P(x)  -  1  -  EXP  t-j  /  Y 

¥  -00 


X  +  X, 


d  x,] 


Note  that  P(x)  Is  by  no  means  a  density  function. 


EXAMPLE  I 


LATERAL  RANGE  CURVE  FOR  THE  DEFINITE  RANGE  LAW: 


In  this  case: 


Y(r) 


^<r0 


0  otherwise 


recalling  from  Figure  1 1 -5 


Therefore : 


F(x) 


*  1  f  x<rQ 


0  If  x>r. 
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And 


1  If  x<  r0 

P(x) 

0  If  x>  rQ 

and  Its  graph: 


Figure  I I -5 

EXAMPLE  II 

LATERAL  RANGE  CURVE  FOR  THE  INVERSE  CUBE  LAW. 
The  detecting  rate  Is 

Y<r)  •  4 

r4 
or 


P(x)  »  1  -  EXP  [- 

Vx* 
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and  Its  graph 


Vt<v^vs 


»  X 

Figure  II-6 

Observe  that  In  this  case  P ( x)  depends  on  the  speed,  and 
assuming  all  other  factors  remain  constant,  the  probability 
of  detection  at  any  lateral  range  decreases  as  the  speed 
Increases. 

As  It  has  been  seen  before  In  two  example.  It  Is  not  very 
difficult  to  obtain  a  compact  mathematical  form  for  P(x). 
Furthermore  the  parameters  Involved  ("r0"  In  example  I  and 
"k"  In  example  II)  can  be  estimated  by  performing  the  experi¬ 
ment  described  In  chapter  I.  Now,  In  a  search  situation,  the 
target  lateral  range  Is  not  known  for  sure.  For,  if  the 
opposite  Is  true,  then  no  longer  Is  necessary  to  search.  To 
take  this  uncertilnty  In  account,  a  random  variable  X,  Is  used 
to  described  th,»  lateral  range. 

Suppose  that  It  Is  desired  to  compute  the  probability, 

P,  of  detecting  a  target  In  a  given  search  no  matter  how  long 
It  takes.  Recall  that  p(x)  was  defined  as  the  probability  of 
detecting  a  target,  given  that  its  lateral  range  Is  x,  and 


observe  that  the  detection  can  occur  at  any  point  on  the 
Infinite  straight  line  path  used  to  define  X,  and  thus,  by 
any  time.  Further,  if  X  Is  a  random  variable,  and  its  density 
Is  denoted  by  fx(x)  ,  then,  from  probability  theory: 

00 

P  -  /  P{x)  f  (x)  dx  (II-XIV) 

0  * 

Which  Is  a  very  useful  expression,  as  the  next  examples  will 
show. 

EXAMPLE  3 

Assume  that  X  Is  uniformly  distributed  In  the  Interval  (0,0). 

This  Is  the  case  of  a  barrier  with  length  2D,  in  the  center 

of  which  the  sensor  Is  located.  If  a  target  tries  to  trespass 

such  a  barrier,  and  Is  equally  likely  to  cross  It  at  any  point, 

then  the  distribution  of  X  Is  the  one  mentioned.  Or: 

r 

i  If  0<x<D 

fx(x)  <  0 

j  o  otherwise 

Where  "D"  stands  for  a  distance  such  that: 

P(x)  *  0  If  x>D  (A) 

what  means  that  the  detection  system  Is  not  effective  at 

ranges  larger  than  0.  In  this  case  (II-XIV)  Is: 

P  *  l  P(x)  dx 

or  what  amount  to  the  same  due  to  equation  (A): 

P  “  i  /  P(x)  dx 
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But,  the  Integral  In  the  formula  above  Is  just  half  of  the 
area  below  the  lateral  range  curve.  By  definition,  this  area 

Is  called  "sweep  width",  or  In  symbols 

«•  00 

W  -  2  /  P (x)  dx  *  /  P(x)  dx 

o  -« 

so  the  formula  for  the  probability  of  detection  becomes 


Observe  that  the  formula  above  justifies  the  use  that  way 
made  of  "w"  In  I.  For,  If  2D  Is  distance  such  that  P*1  the 
definition  of  "w"  used  In  the  mentioned  example  Is  reaflrmed. 
Figure  7  below  pictures  the  situation. 


& 
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/ 
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Figure  1 1 -7 

Observe  that  the  dashed  area  Is  swept  with  efficiency 
w/2D  and  that  Is  the  motive  of  the  name  sweep  width.  Under 
the  same  assumption  about  the  lateral  range  distribution, 
two  expressions  for  "w",  will  be  derived: 


SWEEP  WIDTH  FOR  DEFINITE  RANGE  LAW 


In  this  case,  recalling  from  example  1: 
P(x) 


/- 

(o  ot 


x£ro 


otherwl se 


Then 


W  -  2 


2  ro 


t 


I 
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V 


which  Is  the  sweep  width  used  as  example  In  I.  It  Is  also 
twice  the  dashed  area  In  Figure  5. 

EXAMPLE  5 

SWEEP  WIDTH  FOR  INVERSE  CURVE  LAW 
In  this  case  as  established  In  example  2: 

P(x)  -  1  -  EXP  [- 

Vx‘ 


Then 

oo 

W  -  2  /  P(x)  dx  « 

0 

oo 

*  2  /  {1  -  EXP  [-  -?£]}  dx 

o  Vx‘ 

The  result  can  be  obtained  by  the  use  of  Integration  by  parts 
W  »  2  /TW 

nr 

In  real  life  situations,  the  track's  lengths  are  not 
infinitely  long.  However,  If  they  are  long  enough  In  the 
sense  that  the  sensor's  detection  capability  tends  to  zero 
beyond  a  certain  point  In  the  path,  then  the  preceding  dis¬ 
cussion  applies.  Suppose  the  length  of  the  track  Is: 

Y  -  xj  -  x' 
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such  that 


y(r)  -  0 


If  r  >  /y 


x£  -Tv1  or  r  *  '777^7 

For  any  lateral  range  x,.  For  this  kind  of  situation  the  time 
during  what  the  target  Is  exposed  to  detection  Is: 


T  « 


Because  after  this  time  yCr ( t ) 3  *  0»  then,  the  detection  rate 
along  the  path: 

x*  f 


but 


P  -  J- 
y 


so 


’  IDT 


or,  substituting  the  value  of  "T" 


i  -  MV 

x  ?57 


but  2DY  Is  the  area  swept  by  the  sensor,  then 


i  b  wv 

X  B  A” 


(II-XV) 
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Is  the  average  detection  rate  along  the  considered  path.  Now, 
suppose  a  sensor  with  velocity  and  sweep  width  known  and  a 
target  whose  location  has  a  uniform  distribution  In  the  search 
area.  If  It  Is  known  how  many  paths  are  Intended  to  be  placed 
In  the  S.A.  and  the  time  spent  In  each  one  of  them,  then  It  Is 
possible  to  compute  the  probability  of  finding  the  target  by 
a  suitable  combination  of  each  track's  probability,  which  Is 
a  function  of  X.  Considering  that  the  sum  of  the  time  spent 
In  each  path  can  be  easily  transformed  In  a  total  length  and 
then  In  a  density  of  effort,  that  probability  of  detection 
can  be  written  In  terms  of  density  of  effort,  which  Is  just 
the  definition  for  L.D.F.  Consider  the  following  examples: 

EXAMPLE  6 

L.D.F.  FOR  RANDOM  SEARCH 

Suppose  that: 

(a)  the  target's  location  Is  uniformly  distributed  Inside 
an  area  "A"  which  Is  extremely  large  compared  to  the  sensor's 
Influence  area. 

(b)  The  sensor's  vehicle  movement  Is  such  that  any 
region  of  the  area  Is  equally  likely  to  be  covered. 

(c)  The  tracks  are  mutually  Independent. 

(d)  The  search  vehicle's  speed  Is  always  the  same. 

The  preceding  assumptions  characterizes  a  random  search. 

For  this  case  the  detection  rate  Is  a  constant  at  any  Instant 
In  time.  For,  due  to  the  Independence  of  the  paths,  the 
Instantaneous  probability  of  detection  Is  the  same  for  them  all. 
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By  a  reasoning,  which  Is  similar  In  all  aspects  to  the 
one  presented  In  II-E  for  the  sighting  potential's  cases,  the 
probability  of  detection  at  or  before  time  "t"  Is 

P(t)  -  1  -  EXP  [-Xt]  -  1  -  EXP  [-^] 
which  leads  to  the  L.D.F. 

B(z)  «  1  -  EXP  [-WZ]  (II-XVI) 

that  was  already  derived  In  I. 

EXAPLE  7 

L.D.F.  FOR  SYSTEMATIC  SEARCH 

The  general  case  of  systematic  searching  by  parallel 
sweeps  Is  treated  In  [1].  Now,  here  will  be  derived  the  case 
when  "W"  Is  given  by  a  definite  range  law  and  stated,  without 
proof  (which  also  can  be  found  In  [1]),  the  result  for  the 
Inverse  cube  law.  For  the  first  case  recall  that  W*2rQ  and 
that  within  such  a  distance  the  detection  Is  sure  to  happen. 
If  “A"  stands  for  the  total  search  area,  the  portion  of  "A" 
which  was  not  yet  swept  after  time  Ht"  Is 

A  -  WVt  If  0<t<A 

A(t)  •  "W 

0  If  yy  <t<* 

and  by  this  time  the  detection  rate  Is 

»<*>  ■  STTT 
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therefore  the  probability  of  detection  by  time  "t"  Is: 

1  -  EXP  [-  /  A ( t ) d t ]  *  If  o<t<gy 
o 

1  ^f  <t 

which  leads  to  the  L.O.F. 

If  o<t<£y 

If  ay  <t 

that  agrees  again  with  the  one  derived  In  I.  For  the  second 
case,  the  Inverse  cube  law,  the  sweep  width  seen  In  example  4 
was  found  to  be 

W  *  2  /2t rk 

T“ 

Through  some  more  elaborate  arguments  than  the  preceding  ones, 
but  still  relying  In  the  same  general  principles,  It  can  be 
shown  that  the  L.D.F.  for  the  Inverse  cube  law  Is 

B(Z  )  *  2  ♦  (jp-  WZ)  -1  (II-XVIII) 

Where  4>  stands  for  the  cumulative  distribution  function  of 
the  standard  normal  di  strl butl on .  Now  It  Is  possible  to  draw 
Figure  6,  chapter  I  again.  Including  the  L.D.F.  (II-XVIII). 
This  Is  done  In  Figure  II-8  below. 
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Figure  1 1 -8 

As  was  said  before  In  I-A,  the  best  situation  Is  systematic 
search  with  definite  range  law,  and  the  worst  situation  is  the 
random  search  case.  The  actual  L.D.F.,  when  performing  a 
systematic  search,  must  be  situated  between  them. 

F.  THE  DISTRIBUTION  OF  THE  TARGET'S  LOCATION 

The  results  established  In  the  anterior  sections  depend 
In  the  uniform  distribution  of  the  target's  location. 

However,  It  Is  not  always  the  case  that  this  type  of  distribution 
can  be  assumed  to  describe  the  target's  position.  The  way  to 
overcome  this  difficulty  Is  by  means  of  a  search  grid  which  Is 
constructed  as  follows: 

Suppose  a  bi dlmensi onal  distribution  Is  adopted  to  describe 
the  target's  location  and  that  its  C.D.F.  Is  denoted  by  Ffx-j.Xg)* 

[x, »x2] 
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Then  the  probability  that  the  target  Is  located  In  cell 
Rk  of  Figure  II-9  Is 

Pk  *  F^.Xg)  -  F(x1  .Xg)  -  F(x^  ,x!|)  +  Ffx^  .xj  )  (II-IXX) 

[x^]  Cx,  x23  [x,  x2]  [xltx2] 


Figure  1 1 -9 

Divide  the  area"A"  under  consideration  Into  "J"  nonToverlapp1  ng 
cell  In  such  a  way  that: 

J 

l  P,  ■  1 

M  J 

If  we  assume  that  for  any  "j"  the  probability  "P j "  Is 
uniformly  distributed  Inside  the  cell  "Rj",  then  for  each 
cell  all  the  preceding  discussions  apply.  The  mentioned 


assumption  leads  to  replace  the  original  continuous  distribu¬ 
tion  by  a  discrete  one.  If  the  size  of  the  cells  Is  not  too 
big,  the  new  model  Is  not  very  different  from  the  actual  one. 

It  Is  a  very  common  practice  to  center  the  grid's  coordinate 
system  In  the  highest  probability  cell  and  set  the  axis  parallel 
to  the  north/south  orientation  of  the  navigation  chart.  If  the 
domain  of  the  bi dlmensl onal  distribution  Is  Infinite,  the 
S.A.  encloses  the  desired  probability  to  be  divided  In  the  cells. 
Evidently  In  this  case,  the  cells'  probabilities  do  not  quite 
add  up  to  1.  Further  details  are  found  In  [7]  and  [8].  Some 
examples  of  distributions  most  often  encountered  In  search 
problems  will  be  given  below: 


EXAMPLE  7 

BI DIMENS  I ONAL  UNIFORM  DISTRIBUTION  (INDEPENDENT  X1  AND  X2) 
Although  this  Is  a  very  well  known  distribution.  It  is 
worth-while  to  state  Its  formal  definition: 

- . v - lb 


Consider  the  rectangle  [a,b,c,d]  of  Figure  10.  A  vector 
random  variable  [X^ ,  Xg]  has  a  bldlmenslonal  uniform  distribu¬ 
tion  defined  over  [a,b,c,d,]  if  its  density  function  Is: 


f  (xltx2) 
CX1 ,X2] 


The  C.D.F.  Is 


(*1  -  xl>  <x2  *  *2) 


If  *1  "  X1  *  x) 


X2  <  X2  <  X2 


0  otherwise 


Xi  X  p 

F(x-|  ,x2)  «  /,  /,  f  (xltx2)  dx-j 

tW  x’  Xl  [*,.*.] 

or  be  performing  /.he  intergration  1  c 


‘2 


F ( x 1  »x2 )  = 

[X^g] 


0  1j  *,  <  x,  or  x2  <  x2 


(  X^j  “X]  )  (x2-x2) 

- n - 1 - n - 1 — 

( x •j  -  x})  (x2  -  x2) 


x1  <  <  x] 


x2  <  x2  <  x2 


^  It 

/X1 >  X1 

P/  and 

\x2  >_x" 


Observe  that  In  Figure  10,  the  C.D.F.  argument  Is  repre¬ 
sented  by  the  dashed  region  and  its  value  is  given  by  the 
ratio  between  the  region's  area  ;nd  the  total  area,  as  the 
name  of  "uniform"  intuitively  suggests.  Another  point  is  that 
the  marginal  densities  for  X1 ,  and  X2  are: 


Py  (x,)  *  /  *  f  (x,M)  du  = 

X1  1  xA 


[x1,x2] 


(x1  -  x]  ) 
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0  • 


f-  (x,)  ■  Z,1  f  (H  X2)  du  -  - - - 

1  X1  [X1 Ix23  (X2  ‘  x2> 

If  x^  and  x2  are  Inside  the  proper  range  and  "0"  otherwise. 


Then 


f  (x1 ,x2)  «  f  (x1 )  fx  (x2) 

rw  w  .  X1  2 


[X,,x2] 


and  also 


F  (x,,x2)  *  F^tx,)  FX2  (x2) 

[Xj.Xj] 

all  considering  the  proper  range,  r  course.  The  above  results 
characterize  the  Independence  of  the  random  variables  X-|  and 


EXAMPLE  8 


THE  BIVARIATE  NORMAL  DISTRIBUTION 


Th4s  is  a  very  useful  distribution  In  real  life  problems, 
for  It  Is  a  well  known  fact  that  error  distributions  are  In 
general  normal.  Navigation  error  can  be  decomposed  In  two 
components:  Latitude  and  longitudo,  and  If  both  are  univariate 
normal  their  composition  Is  bivariate  normal.  Example:  If  a 
ship  In  distress  asks  for  help  and  gives  her  last  fix,  It  Is 


reasonable  to  use  this  distribution  to  describe  her  position. 

By  using  x-j  to  denote  the  east/west  direction,  Xg  for  the 

north/south.  M  for  mean  6  for  standard  deviation  and  pY 

a2 

for  coeflcient  of  correlation,  the  bivariate  normal  density: 
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f  (x1  ,x2) 


1 


[xrx2] 


2iroy  Ow  (1-p  ) 
*1  *2 


EXP  {- 


[(■ 


xl"Mx 


2(l-pTf '  °X1 


i)2 


2p  (■ 


xrMx, 


-)  (■ 


VMx, 


-)  +  (■ 


x2-Mx 


i)2]} 


If  however  It  Is  assumed  that  the  north/south  and  east/west 

errors  are  independent  (which  is  common  practice)  then  p=0. 

Furthermore  if  the  distribution  Is  centered  in  the  datum  then 

MY  *  MY  *  0.  With  these  simplifications  the  above  formula 
*1  *2 
becomes: 


f  (x1  ,x2) 

[xrx2] 


1 

2lT0y  Oy 

*1  *2 


EXP  {- 


1  r 

*1 


Figure  11  shows,  the  shape  of  f(XjX2)  and  the  contours  of 

[X^] 

equal  probability,  which  are  generally  ellipses  with  the  main 

axis  laying  In  the  coordinate  axis,  which  corresponds  to  the 

higher  variance.  The  contours  are  circles  when  6W  *  6V 

X1  x2 


Considering  that  the  values  of  the  standard  normal  CDF 
are  easily  available  In  many  tables  and  that,  due  to  their 
Independence 

F  (*,,x2)  -  F^  (x,)  Fj^(x2) 

[*,.x2] 


It  Is  very  simple  to  compute  any  region's  probability, 
particular  case.  In  reference  to  Figure  I I -9 ,  formula 
becomes: 


[F 


X^xl^  " 


Vxi)] 


[Fx2(x2>  ‘  FX2(x2^ 


In  this 
(II- 1 XX) 


But 

F«(x)  -  *  <**) 

In  this  case  M*0,  then 
Fx(x)  -  *  (S-) 
and 


P 


K 


c*  (^-)  -  ♦  (^-)]  [*  < 

aX,  aX, 


II 


*  (=H3 


(*) 


Where  *,  Is  the  standard  normal  CDF.  In  order  to  derive  the 
dimensions  of  a  search  area  for  this  case,  the  common  usage  is 
to  construct  a  rectangle,  centered  on  the  datum,  whose  sides' 
length  are  usually  expressed  in  terms  of  standard  deviation, 
and  that  encloses  the  desired  amount  of  probability.  This 
can  be  done  by  replacing  In  (*)  P^  by  the  desired  probability, 
and  solving  the  equation  backwards,  taking  advantage  of  the 
symmetry  of  the  normal  distribution.  Figure  12  picture  the 
case  of  a  search  area  for  a  circular  normal  distribution,  In 
this  case  (*)  becomes: 
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I  III 

Xi  x«  x«  x« 

Pk  -  [♦  (-£)  -  4  (-  C*  -  *  (  o^  * 
x!  '  2 

C*  (4>  '  ♦ 


Figure  11-12 

And  by  consulting  a  table,  the  solution  Is  straight  forward. 
As  Illustration,  some  values  are  tabled  below: 


Side  Length  Probability  Enclosed 

1  Standard  Deviation  0.466 

2  Standard  Deviation  0.911 

3  Standard  Deviation  0.994 

4  Standard  Deviation  “  1.0 


Finally  there  is  another  way  of  writing  the  density  of 
this  distribution  which  sometimes  Is  very  helpfull.  It  Is 
the  polar  from  and  the  formula  that  will  be  stated  holds 


for  the  circular  case  (6V 

X1 


„  r2 

f  (r ,0)  ■  EXP  [-  -V) 

2tt  ac  2cr 

C  R-|  ©3 
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Where,  as  usual 


Xj  ■  r  cos  6 

x2  ■  r  sen  0 

r  ^  0 
0  <  0  <  2tt 

EXAMPLE  9 

MULTIPLE  SCENARIOS  DISTRIBUTIONS 
An  Interesting  problem  which  arises  In  searches  on 
large  scale,  (for  example  the  search  for  the  H-bomb  dropped 
accidentally  In  the  Spanish  cost  In  January  1966,  see  [7] 
for  details)  Is  the  problem  of  multiple  scenarios. 

Suppose  that  from  two  different  sources  of  Information, 
two  different  datums  are  given  to  the  analyst  In  charge  of 
planning  the  search.  Figure  13  pictures  the  situation. 


,  •  / 
/ 


Figure  11-13 


i 


Even  assuming  a  circular  normal  distribution*  the  problem 
of  where  to  center  it  still  remains.  The  best  approach,  Is 
to  center  a  distribution  on  each  datum  and  by  weighting  the 
reliability  of  each  source  by  a  number  between  0  and  1,  CRJ 
(which  reads  "credence  In  source  "J"),  obtain  a  distribution 
which  then  Is  used  as  the  target's  location  distribution. 
Consider  Figure  13  and  suppose  that  the  probability  In  the 
dashed  cell  Rk  Is  to  be  computed.  From  probability  theory 
Is  recalled  that: 

Prob  (Rk)  =  Prob  ( Rk/DATUM  I  Is  CORRECT)  Prob  (DATUM  I 
is  CORRECT)  +  Prob  ( Rk/DATUM  II  Is  CORRECT)  Prob  (DATUM  II 
Is  CORRECT) 

IF  CRJ  Is  taken  as  the  probability  of  datum  j  being  correct, 
then: 

P(Rk)  -  P(Rk/I)  CRI  +  P(Rk/II)  CRI I 
Care  should  be  taken  about  the  value  of  CRJ.  This  Is  a 
very  subjective  matter,  and  can  have  an  enormous  Influence 
In  the  search. 

EXAMPLE  10 

OTHER  DISTRIBUTIONS 

Although  the  normal  and  the  uniform  are  the  most  usual 
distributions  In  search  problems,  there  are,  besides  the 
multiple  scenarios  situations,  other  cases  where  they  do 
not  apply,  at  least  In  the  straightforward  way  presented 
here. 
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Suppose,  for  Instance,  that  a  position  Is  available  from 
a  target,  say  normally  distributed  around  a  datum,  and  a  radio 
direction  finder  bearing,  with  a  uniformly  distributed  error 
of  k  degrees  to  each  side.  The  situation  can  be  visualized 


It  Is  not  hard  to  Imagine  many  other  analogous  situations. 
In  any  case,  however,  the  procedure  Is  the  same,  combining 
the  probabilities  by  conditioning,  to  get  a  search  grid. 
Sometimes,  a  computer  Is  needed,  either  because  of  the  big 
volume  of  operations  or  even  due  to  the  analytical  Impossi¬ 
bility  of  deriving  the  conditional  distribution.  In  this 
last  case,  a  monte  carlo  method  must  be  used  (see[7]). 

\  * 


G.  SEARCH  PLANS 

In  I  the  Ideas  of  allocation,  and  of  the  cost  of  an 
allocation  were  discussed.  Here,  a  more  general  Idea  of  a 
search  plan  Is  Introduced.  Consider  a  function  M(j,<f>(t)) 
such  that: 

(a)  M(j,<t>)  Is  a  allocation  for  each  $ 

(b)  M(j,$)  Is  Increasing  In  <J>  for  j*l....J. 

This  function  Is  called  a  search  plan.  The  replacement  of 
j  by  a  vector  (xj>  x2),  leads  to  the  definition  for  the 
continuous  case.  Condition  a)  means  that,  for  each  cell 
(point)  the  function  H  assigns  a  density  of  effort  z,  which 
depends  also  In  4  ,  the  amount  of  effort  available.  Condition 
b)  is  the  formal  statement  of  the  fact  that  It  is  Impossible 
to  take  effort  out  of  a  cell  (point),  once  it  was  placed  there. 

The  definition  of  cost  of  an  allocation  carries  over  to 
the  case  of  search  plans.  In  fact,  for  the  discrete  case  It 
Is  the  same  as  defined  In  I,  while  for  the  continuous, 
the  sum  Is  replaced  by  an  Integral,  an  It  becomes: 

C(M(x,4>)  *  /  M(x,<J>)dx 

“  x 

Suppose  now  that  <J>  Is  an  Increasing  function  of  the  time, 
♦=4>(t).  The  problem  which  was  solved  In  I,  with  the  track 
length  measuring  the  effort,  and  given  by 

Vt  =  L 


where  V  remained  constant,  was  an  example  of  this  situation. 


In  this  case,  $(t)  Is  simply  Vt.  If  the  search  plans  Is  such 
that 

a)  C(M( j ,t) )  -  *(t) 

b)  Prob  (M( j ,t ) )  -  MAX  {Prob  (2) :C(Z)  <  *  (t)} 

Then  M  Is  a  UNIFORMLY  OPTIMAL  SEARCH  PLAN,  denoted  M#. 

Observe  that  the  above  definition  formalizes  what  was 
Intuitively  shown  In  I.  For,  If  at  any  time  t,  the  plan 
costs  exactly  the  amount  of  effort  by  that  time  available, 
and  reaches  the  maximum  value  of  probability  of  detection 
which  Is  attainable  with  that  effort,  then,  besides  of  being 
the  best  possible  plan  for  the  constraint  which  Is  Imposed, 

It  Is  "constraint  Independent"  In  the  sense  that  no  constraint 
alteration  will  be  able  to  alter  the  uniform  optimum  search 
plan. 

In  the  next  section  It  will  be  discussed  how  to  construct 
uniformly  optimal  search  plans. 

H.  RETURN  FUNCTION  AND  UNIFORMLY  OPTIMAL  SEARCH  PLANS 

According  to  what  was  stated  In  I,  the  mathematical 

formulation  of  the  resource  allocation  problem  Is: 

MAX  P (2 ) 
s.t. 

C  ( 2 )  <  C 

Which,  In  the  context  of  uniformly  optimal  search  plans,  is 
equivalent  to: 

MAX  P  (M(j,*)) 
s.t. 

C  (M(j,*))  <  *  (t) 
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Therefore,  the  problem  is  to  find  a  search  plan  which  by 
spending  at  each  Instant  Mt"  all  the  effort  which  by  that 
time  Is  available,  maximizes  the  probability  of  detection. 

The  discussion  of  the  general  solution  for  this  type  of 
problem  depends  strongly  on  mathematical  programing  methods, 
and  Is  far  beyond  the  scope  of  this  paper.  The  Interested 
reader  can  see  [2],  and  [9],  for  applied  and  theoretical, 
respectively,  details.  However,  It  is  possible  to  devise  a 
solution  for  the  mentioned  problem,  by  using  the  concept  of 
a  return  function,  whose  meaning  can  be  Intuitively  explained, 
as  following: 

Suppose  that  an  Increment  of  effort,  h.  Is  applied  In  cell 
j.  The  corresponding  Increment  In  probability  of  detection, 
considering  that  an  amount  z  was  already  applied.  Is: 

Pj  [ B j  (Z  +  h)  -  (Z)] 

and  the  Increment  In  cost  Is 

hA. 

but 

Pj  [Bj  (Z  +  h)  -  (Z)]  :  Pj  b]  (Z)  h 

therefore 

s 

Is  the  ratio  between  the  Increase  In  effectiveness  and  the 
Increase  In  cost,  that  results  from  placing  a  small  Increment 
of  effort  In  cell  j.  This  ratio  is  called  the  return  function, 
and  Is  denoted  by  A(z).  But  it  Is  also  a  H.O.E.  of  what  Is 
being  done,  for  the  mentioned  Increment  of  effort  will  produce 


Aj h 


better  results  where  x(z)  Is  maximum.  In  fact,  (see  [2]  for 
proof)  It  can  be  formally  stated:  If  M  Is  the  search  plan 
such  that  the  next  Increment  of  effort  Is  applied  In  the  cell 
where  X(z)  Is  highest,  then  M=M*1s  the  uniformly  optimal  plan. 
The  preceding  discussion  which  has  been  done  for  the  discrete 
case,  applies,  with  suitable  transformations,  to  the  continuous, 
as  It  will  be  shown  In  the  examples, 

EXAMPLE  11  (discrete  case) 

RETURN  FUNCTION  FOR  THE  RANDOM  SEARCH  L.D.F. 

Let 

Bj  (Z)  =  1  -  EXP  [-  WZ]  ,  j  -  1,  2  ...  J 

then 

Bj  (Z)  *  EXP  [-WZ] 

and 

P,  EXP  [-WZ) 

\i  (z)  *  - 

Suppose  that  all  the  cells  have  the  same  area.  Then, 
the  next  Increment  of  effort  must  be  placed  In  the  cell  where 
Xj(z)  Is  highest.  But  this  Is  exactly  the  procedure  used  In 
the  heading  I  example,  for  Pj  EXP[-wz]  Is  the  numerator  of 
the  posterior  probability  formula  for  cell  j,  and  due  to  the 
equality  of  the  denominators,  nly  the  numerators  were  taken 
In  to  considerations. 

Therefore,  for  the  considered  L.D.F.  both  procedures 
are  equivalent.  Furthermore,  this  equivalence  Is  a  proof  that 
the  search  plan  used  In  the  mentioned  example,  was  uniformly 
optimal . 


76 


EXAMPLE  12 

AN  UNIFORMLY  OPTIMAL  PLAN  FOR  A  CONTINUOUS  SEARCH  AREA  CASE 


Suppose  the  situation  described  by  the  following  assump 
tlons : 

a)  The  target's  location  Is  described  by  a  circular 
normal  distribution,  1  ,e: 


f  (r.8) 
[R-l  0] 


r  ^  0 

0  <  0  <  2tt 

b)  The  LDF  for  the  case,  which  Independs  In  (Xj,  x?) ,  Is: 


B  (Z)  =  1  -  EXP[-WZ] 
and,  therefore,  assuming  W*1 
B * (Z)  *  EXP  [-Z] 


c)  The  available  search  effort  Is  given  by 


♦(t)  «  WVt 


d)  The  return  function,  for  the  continuous  case.  Is 
defined  as: 


A(x,Z)  =  f  (x1  ,x2)  B  (x,Z) 


Cx,  ,x2] 

or,  In  polar  coordinates: 

A(r,Z)  =  f  (r,0)  B  ( ( r] 0  ,Z) 

[R.e] 
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for  a  given  z. 


To  obtain  an  uniform  optimal  plan.  It  Is  necessary  to 
find  X  for  each  time  t,  namely  X(t)  and  then  solve  this 

equation  for 

Z  -  M1  ((r.e).t) 

because  M*((r9),t)  )  Is  the  amount  of  effort  that  must  be 

placed  at  point  (r^0),  by  the  time  t. 

Under  these  assumptions.  It  Is  possible  to  derive  an 
expression  for  X ( t )  (see  [2]  for  details)  and  the  expression  Is 

x(t>  •  —V  exp  t- 

2tto  no 


by  letting 


the  search  plan  Is: 


where 


( (r1e)1t)  =  [In  (f  (r1e)/X(t))] 


.  [K/t  -  -^]+ 
Zac 


[X]+  * 


X  If  X>0 
0  If  X<0 
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which  does  not  depend  on  d.  This  Is  due  to  the  fact  that 


f(t*i  )  Is  radially  symmetric  and  also  because  B(Z)  does  not 

C^e] 

depend  on  position.  Then*  supressing  0,  It  Is  possible  to 
write  the  plan  In  a  more  convemient  way: 


K/t  - 


r2 

i? 


For  r  <  R(t) 


0 


For  r  >  R(t) 


where 

R( t)  *  2o2  K/t  ,  t  >  0 

Therefore,  the  optimal  plan  starts  at  the  origin  and 
spreads  in  such  a  way  that  by  time  t  all  the  effort  Is  placed 
Inside  a  circle  of  radius  R(t).  This  means  that  as  the  search 
progresses,  the  effort  Is  accumulated  around  the  origin. 

This  makes  sense  because  the  origin  has  the  highest  prior 
target  location  density. 

Such  a  plan  can  only  be  approximated  in  real  life  situations. 


I.  EXPECTED  TIME  FOR  FINDING  THE  TARGET 

In  I,  It  was  pointed  out  that  the  time  taken  to  find  a 
target,  can  be  used  as  a  M.O.E.  for  search  plans.  However, 
it  should  noticed  that  this  time  is  a  random  variable  T. 

A  reasonable  approach  Is  to  minimize  the  expected  value 
of  the  random  variable  T.  Let  y  be  the  expectation  of  T,  and 
let  P(t)  be  the  probability  of  finding  the  target  at,  or  before 
time  t. 
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Notice  that  P(t)  Is  the  C.D.F.  of  T.  By  definition,  the  value 
of  y  Is : 

p  *  /  t  P  (t)  dt 
-00 

but,  due  to  the  fact  that  T  Is  a  non  negative  random  variable, 
It  can  be  demonstrated  (see  [10]),  that: 

V  I  [1  -  P(t)]  dt 
o 

Now,  consider  the  uniform  optimal  search  plan  M*.  By  the 

§ 

very  definition  of  M  : 

P[M#(x,t)]  >  P[M(x,t)] 

where  M  Is  any  other  search  plan.  Therefore: 

1  -  P[M#(xft)]  <  1  -  P[M(x ,t)3 
and 

0 9  M  90 

f  (1  -  p  [M#(x,t)]}  dt  <  /  (l-P[M(x,t)]}  dt 

0  0 

or 

<.  M 

§ 

where  y  and  y  are  the  mean  times  to  find  the  target,  for 

I 

M  and  M,  respectively. 

Therefore,  the  uniformly  optimal  search  plan  minimizes 
the  expectation  of  T,  and  under  these  circumstances  it  is 
imaterlal  what  M . 0 . Z . 


is  used  In  the  problem. 


EXAMPLE  13 

EXPECTED  TIME  FOR  THE  UNIFORM  TARGET'S  LOCATION  DISTRIBUTION 
Suppose  that  the  following  conditions  holds: 

(a)  The  target's  location  Is  uniformly  distributed 
Inside  an  area  A. 

(b)  The  LDF  Is: 

B(Z)  -  1  -  EXP  [1WZ] 
therefore,  the  CDF  for  T  Is: 

FT(t)  •  P(T<t)  -  P(t)  *  1  -  EXP  [-Sjp  ] 


under  these  conditions,  the  expected  time  for  finding  the 
target  Is: 

m  m 

E[T]  -  u  -  !  [1  -  P(t)[dt  ■  /  EXP  [-H]£jdt  -  L- 

0  o  M  wv 


Applying  this  formula  to  chapter's  I  example.  In  the  case 
Pj  ■  1/4  for  all  cells,  and  recalling  that  the  total  area 
has  6.000  square  miles,  the  expected  time  Is: 


y 


6.000 

15x100 


4  hours 


Observe  that  If  In  the  CDF,  t  Is  replaced  by  y,  the  probability 
of  detection  after  the  expected  time  Is  gone  Is: 

UVA 

P(y )  -  1  -  EXP  [-{fi/j]  »  1  -  EXP  [-1]  »  0.633 
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The  preceding  computations  were  done  under  the  assumption 
that  the  search  only  stops  If  the  target  Is  located.  However, 
It  Is  possible  to  compute  the  expected  time  to  detect  the 
target,  under  the  assumption  that  the  search  continues  either 
until  the  target  Is  located,  or  until  some  pre-asslgned  value 
of  probability  of  detection,  say  P*,  Is  reached.  Let  T* 
be  the  time  when: 

P  (t)  -  P* 

from  condition  b) 

P  (T*)  -  P*  ■  1  -  EXP  ] 

and  therefore 

T* '  w  1n 

Let  e be  a  random  variable  which  describes  the  duration  of 
such  a  search.  The  C.D.F,  for  6,  can  be  derived  considering 
that: 

P  (0<0  if  0<T*)  -  1  -  EXP 


and 


P  (0<0  if  0>T*)  -  1 


Then 


f  (e) 


1  -  £XP  [~  If  0<T* 

l  if  e>T* 


and  the  expected  value  of  9  Is: 


•  T*  UUQ 

E  [e]  »  /  (l  -  f  (e) )  de  -  /  exp  [-  -*-]  de 

O  0 

♦  /  (l-i)  de  -  /  exp  [-  MV©]  de  • 

T*  0  A 


Suppose  that  In  the  6.000  square  miles  area  of  Heading' 
I  example,  a  search  Is  Intended  to  be  performed,  which  Is 
to  be  stopped  when  a  probability  of  detection  P*  Is  reached. 
In  order  to  gain  some  Insight  In  the  order  of  magnitudes 
Involved,  T*  and  E[e]  will  be  computed  for  several  values 
of  P*,  as  It  Is  shown  In  the  table  below: 
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p* 

T*  (Hours) 

E £ 8 3  (Hour: 

0.500 

2.772 

2.000 

0.600 

3.665 

2.400 

0.700 

4.811 

2.800 

0.800 

6.437 

3.200 

0.900 

9.210 

3.600 

0.910 

9.624 

3.640 

0.920 

10.100 

3.680 

0.930 

10.612 

3.720 

0.940 

11.230 

3.760 

0.950 

11.980 

3.800 

0.960 

12.872 

3.840 

0.970 

14.020 

3.880 

0.980 

15.648 

3.920 

0.990 

18.420 

3.960 

0.992 

19.312 

3.968 

0.994 

20.444 

3.976 

0.996 

22.084 

3.984 

0.998 

24.856 

3.992 

1  .000 

T*-mo 

4.000 

It  must  be  noticed  the  enormous  disagreement  between 
the  time  to  stop  and  the  expected  time  for  the  search,  as 
P*  Is  Increased.  Suppose  that  In  a  mine  sweeping  problem, 
the  mines  are  randomly  distributed  In  the  area  A  and  the 
sensor  characteristics  are,  as  before  W  ■  15  mile  and 
V  ■  100  knots.  In  such  a  context,  P(t)  can  be  Interpreted 
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as  the  fraction  of  the  total  number  of  mines,  which  was 
detected  by  time  t.  The  tremendous  decreasing  return's 
effect  can  be  seen  In  the  above  table  (column  1  and  2). 

For  example,  to  Increase  the  fraction  of  mines  detected  from 
50X  to  70*,  the  stopping  time  Is  Increased  In  2.039  hours, 
and  to  Increase  that  fraction  from  99.4*  to  99.8*,  the 
tine  Increases  In  4.412  hours. 

Suppose  now  the  same  area  A,  sweep  width  W,  and  velo¬ 
city  V,  as  In  example  13,  l.e.,  6000  square  miles,  15  miles, 
100  knots,  respectively.  However,  the  probability  distri¬ 
bution  Is  no  longer  uniform  but  given  by 


0.4  If  j  ■  I 

0.3  If  j  -  II 

0.2  If  j  «  III 

0.1  If  j  »  IV 


As  In  the  last  part  of  Heading's  I  example.  Recall  that 
a 'uniformly  optimal  search  plan  was  derived  there,  which  can 
be  written  as: 


M#(j,t) 


Vt 

IF  0  <  t  <_  ©-j 

z 1 1 1  »  ZIV  «  0 


Where 


01  *  0.28770  hours 
02  ■  0.8100  hours 
e3  -  2.07945  hours 
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or  replacing  the  values  for  A y  0„,  W,  V,  and  Pj 


0.4  (1  -  EXP  [It]  0  <  t  <  0.28770 


0.4  (1  -  EXP  [-0.28770  +J]  + 


0.3  (1  -  EXP  [-  £)) 


0.28770  <  t  <  1.0977 


P(t)  -  < 


0.4  (1 
0.3  (1 
0.2  <1  -  EXP[-  £]) 


EXP  [-0.6927+  |])  + 

EXP  [-0.405  +£])  +  1.0977  <  t  1  3.17715 


0.4  (1 
0.3  (1 
0.2  (1 
0.1  (1 


EXP  [-1.3927  +£])  + 
EXP  [-1.105  +{■])  + 
EXP  [-0.7  +  £])  + 
EXP  [-  $•]) 


3.77715  <  t 


As  was  mentioned  before.  In  general: 


U  »  /  [1  -  P(t)]  dt 

o 


In  this  case,  however,  due  to  the  form  of  P(t),  this  expression 
becomes : 


0,  01+09 

y  »  /  1  [1  -  P(t) J  dt  +  /  1  1  [1  -  P(t)]  dt  + 

o  0i 


'1 


01+09+09  * 

/  1  2  3  [l  -  P(t)3  dt  +  / 


®l+®2 


®l+®2t03 


n-p(t)]  dt 


0.2877 


Jl-0.8[l-EXP(-t)]dt  =  0.27342  hours 


1.0977 


jl-0.4[l-EXP(-0.2877-  £]-0. 3[1 -EXP(-  £)]  dt  = 


0.2877 

3.17715 


0.5941216 


jl-0.4[1-EXP(-0.6927-  |)]-0.3[1 -EXP(0.405-  |)]  - 


1.0977 


0.2[1 -EXP ( -  £)]  dt  *  0.5972783  hours 

1  1-0. 4[1-EXP(-1 .3927-  J)]  -0 . 3[ 1 -EXP ( -1  .105  -  J-)] 

1  p 


'3.17715 


-  0. 2[1 -EXP( -0. 7  -  J)]  -  0.1  [1 -EXP ( -  £)] 


0.5270036  hours 


Then 


y  *  I|  +  I2  +  I3  +  I4  s  1.9918235  hours 
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The  computation  above  Is  very  cumbersome,  and.  In  fact, 
for  a  larger  number  of  cells  It  may  be  Impracticable.  In 
real  life  situations  however,  It  Is  possible  to  program  a 
computer  to  do  the  work.  Anyway,  the  knowledge  of  the 
expected  time  to  find  the  target.  Is  very  important  In  search 
planning,  and  this  example  shows  how  to  do  It  for  a  grid 
type  distribution. 

If,  as  It  was  done  In  the  preceding  example,  the  search 

Is  to  be  stopped  If  the  target  Is  not  found  by  time  T*,  when 

P (T* )  ■  P*,  a  pre-asslgned  value,  then,  due  to  the  form  of 

P(t),  becomes  much  more  difficult  from  the  computational  point 

of  view,  to  determine  the  values  of  T*  and  E[e].  The  general 

approach  to  the  problem  Is:  (a)  given  P*,  determine  In 

what  Interval  T*  Is.  Invert  the  corresponding  expression  for 

T* 

P(t)  and  obtain  T*;  (b)  compute  /  [l-P(t)]  dt  ■  E [0 ] 

o 

As  example,  suppose  first  P*  «  0.06  considering  that 
P( 0 . 2877 )  *  0.098;  then  T*  belongs  to  the  first  Interval, 
therefore: 

T*  ■  -  In  [1-  ■  0.16252 

or  In  general: 

A|(  p* 

T*  -  V¥  [In  (1- 

because  In  this  case,  duu  to  the  numbers  which  are  being  used. 


for  any  j. 
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Then 


T* 

E  (© )  •  /  Cl-P(t)]  dt 
o 

or  In  the  example 

0.16252 

E  (e )  -  /  (1-0.4  [1  -  EXP  (-t)])  dt  -  0.157  hours 

o 

However,  If  P*  Is  chosen  to  be  equal  to  0.25,  things  are  not 
so  easy.  Considering  that 

P  (1.097)  •  0.298 

* 

then  T*  Is  located  In  the  second  Interval,  and  the  equation 
to  be  solved  Is: 

P  (x)  «  0.4  [1  -  EXP  (-0.28770  -  £)]  +  0.3  [1  -  EXP(-£)] 
and 

T*  -  0.2877  +  x 

which  only  can  be  solved  by  numerical  methods.  After  the 
value  of  T*  Is  known,  the  expected  time  is: 

0.28770 

E  (9)  »  /  (1-0.4  [1  -  EX P ( - x ) }dx  + 

o 

T* 

/  (1-0.4  [l-EXP(-0.2877  -  £)]  - 

0.2E  J 

-  0.3  [1-EXP  (-  £)]}  dx 
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The  above  calculations,  could  be  done  on  a  computer. 

EXAMPLE  15 

EXPECTED  TIME  FOR  FINDING  A  TARGET  WHOSE  POSITION  HAS  A  CIRCULAR 
NORMAL  DISTRIBUTION 

In  example  12,  It  was  derived  an  M*  plan  for  this  case. 

The  probability  of  detection  at  or  before  time  t  under  this 
plan  Is: 

P(t)  «  1  -  (1  +  K  /l)  EXP  [-  k  /f] 

whose  derivation  can  be  found  In  [2].  The  expected  time  Is 
then: 

/  [1  -  P(t)]dt  -  /  [1  ♦  K  /T]  EXP  [-  k  Vt]dt 

o  o 

2 

•  2I£ 

WV 

Again,  fixing  a  probability  P*,  and  stopping  the  search 
either  If  the  target  Is  found  or  If  a  time  T*  such  that 
P(T*)  ■  P*  Is  reached,  leads  to  another  way  of  defining  the 
expected  time.  In  this  case: 

P*  .  p(T*)  -  1  -  (l  +  k  /T*)  EXP  [-  k  /T*] 

equation  which  can  be  solved  for  T*  by  numerical  methods.  As 
before,  the  expected  time  for  a  search  under  these  conditions 
Is : 
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CO 

E  (0)  ■  /  [1  -  F  (0)]  d0  - 

0 

T* 

»  /  {(1  +  K»/t )  EXP  [-K/t]}  dt  = 

0 

«  T*  +  ^  -  2KEXP(-K/TT) 

K  K  K 

where  6.  F  (0)  and  E(e),  have  the  same  meaning  as  In  example  13. 
J.  FALSE  TARGETS 

Suppose  a  detection  device,  a  radar  for  example,  and 
consider  the  problem  of  setting  it  In  sensitivity  control. 

The  more  sensitive  the  instrument  Is  made  the  less  likely  it 
will  miss  a  target.  Then,  for  any  given  distance  the  instanta¬ 
neous  probability  of  detecting  the  target,  given  it  Is  there, 

Is  an  increasing  function  of  the  radar's  sensitivity.  On  the 
other  hand,  as  the  sensitivity  is  made  higher,  the  probability 
of  a  false  alarm  Is  also  Increased.  For  with  higher  sensitivity 
more  noise  Is  presented  in  the  radar's  scope,  and  since  the 
operator  Is  unable  to  differentiate  between  signal  and  noise, 
he  will  claim  more  frequently  the  presence  of  a  target,  where 
in  fact  what  he  detected  was  noise. 

This  is  one  of  the  several  situations  that  arises  during 
a  search  when  it  is  necessary  to  distinguish  between  the  target 
and  false  alarms.  A  member  of  the  set  composed  of  targets  and 
false  targets  in  called  a  "contact."  Once  there  is  a  contact 
in  the  sensor's  display,  It  takes  some  time  to  identify  it  as 
a  target  or  as  a  false  target,  and  even  then  the  identification 
may  be  incorrect. 
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From  the  above  considerations  ,  It  can  be  concluded  that 
the  presence  of  false  targets  should  affect  the  search  plan. 

In  fact,  under  these  circumstances,  the  search  takes  place  In 
two  phases:  The  first,  called  the  broad  search,  Is  conducted 
using  a  sensor  which  Is  able  to  detect  the  contacts,  but  can 
not  Identify  them.  The  second,  the  contact  Identification 
phase,  Is  conducted  by  another  type  of  sensor  (In  most  cases, 
the  human  eye),  which  Identifies  the  contacts. 

At  this  point,  It  Is  necessary  to  make  the  following 
assumptions  In-order  to  build  a  model  to  handle  the  problem: 

(a)  The  search  area  Is  divided  Into  J  cells,  each  one  with 
area  A j ,  and  having  a  false  target  density  6j.  The  actual 
number  Nj  of  false  targets  In  Aj  Is  Poisson  distributed: 


P(Nj  »  n) 


(Aj^EXPNA^) 


ni 


Let  Tj  stand  for  mean  time  to  Identify  a  contact  In  cell  j. 

(b)  As  soon  as  a  contact  Is  detected,  the  broad  search  stops, 
and  the  contact  Is  Investigated  until  positive  Identification. 

If  It  Is  Identified  as  a  false  target,  the  broad  search  starts 
again. 

Under  the  preceding  assumptions.  It  can  be  demonstrated, 
that  It  Is  not  possible  to  find  a  uniformly  optimal  search  plan. 
However,  It  Is  possible  to  devise  a  plan  which  minimizes  the 
expected  time  to  find  the  target  (observe  that  In  this  case  the 
two  M.O.E.  do  not  coincide). 
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The  following  results  will  be  stated  without  proof 
(complete  discussion  can  be  found  In  [2]): 

(R.  1 )  The  return  function,  for  this  case  Is  given  by: 

A,(Z)  - 

i  Aj[l  ♦  SjTjBjtZ)] 

(R.2)  The  best  search  plan,  M*,  Is  the  one  which  allocates 
the  next  Increment  of  effort  to  the  cell  (or  cells)  where  A(Z) 

Is  higher. 

(R.3)  The  probability  of  finding  the  target  is  still: 

P(Z)  .  j)  Pj  Bj  (Zj) 

However,  when  computing  the  probability  at  or  before  time  T, 

It  Is  to  be  understood  that  this  time  refers  to  broad  search 
time,  not  Including  therefore  the  time  spent  In  contact 
Investigation. 

(R.4)  The  cost  of  an  allocation  Is  In  this  case. 


Where  Zj  refers  to  the  amount  of  broad  search  effort  allocated 
to  cell  j.  For  sake  of  clarification,  suppose  that  Zj  stands 
for  the  time  per  unit  area  spent  at  cell  j,  then: 

J 

C(Z)  =  Z  A.  t,  +  A<  T,  B .  (tj 
j,]  J  J  J  j  J  J  J 

At  this  point,  an  example  seems  to  be  helpful. 
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EXAMPLE  16 

Suppose  that  a  helicopter  Is  searching  for  a  trawler  X. 
whose  location  distribution  Is  pictured  In  Figure  15  below: 


3o  mi 

Somi  So  mi 

Figure  11-15 

Suppose  further  that  there  are  several  other  trawlers  In  both 
areas,  and  the  helicopter  must  approach  each  one  It  contacts 
In  order  to  verify  If  she  Is  the  trawler  X.  If  It  Is  assumed 
that  the  trawlers  are  uniformly  distributed  In  each  area,  the 
model  presented  applies,  playing  X  the  role  of  target,  and  all 
other  trawlers  the  role  of  false  targets. 

Let  the  LDF-be: 

B(Z)  -  Bj  (Ji)  -  1  -  EXP  [- 

where  tj  Is  the  amount  of  time  spent  In  cell  j. 

And  the  other  parameters: 

Tj  =  12  minutes  *  0.2  hours 
V  -  100  knots 
W  »  15  miles 
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And  due  to  past  experience  In  that  fishing  area,  It  Is  estimated 
that  there  are  about  3,000  trawlers  In  cell  I,  and  75  In  cell 
II  then: 

6j  »  2/sq.  ml. 

6j j  ■  0.05/  sq.  ml  . 


In  order  to  know  where  to  start.  It  Is  necessary  to  compute 

the  return  function  for  Z  *  0.  In  this  case: 

WVt, 


■* 


Pj  MV  EXP  (-  -p-) 

3  .  Mvt  — 
AJC1  +  {jTj  MV  ‘XP  <-  -jp)3 


Therefore,  the  plan  starts  In  cell  II,  the  one  with  smaller 
prior  probability.  This  happens  because  of  the  big  difference 
In  false  targets'  density,  over  weighted  the  difference  In 
probability.  The  effort  will  be  placed  only  In  cell  II,  until 
a  time,  0  such  that: 

Rn(9)  =  RjO) 
or 

0.45  EX?  [-0]  _  „  ,„-4 

TV  15  EXP  1-0]  -  9X10- 
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Which  solved  for  gives 


9  *  3.9  hours 

after  this  time,  the  effort  must  be  split  evenly  among  the 
cells. 

K.  SEARCH  PLANNING 

The  following  are  basic  topics  In  search  planning: 

(a)  Determination  of  the  initial  target  location 
probability  distribution. 

(b)  Evaluation  of  the  search  system  capabilities 

(c)  The  decision  to  begin  search 

(d)  The  decision  to  terminate  search 

Before  examining  each  one  of  the  mentioned  topics.  It 
should  be  noticed  that  the  above  list  Is  by  no  means  exhaustive. 
In  fact,  It  Is  very  likely  that  for  each  situation  a  multitude 
of  other  factors  must  be  considered.  However,  as  was  mentioned 
In  the  very  first  sentence  of  this  section,  these  are  basic 
topics  In  the  sense  that  they  are  worthy  of  consideration  under 
any  circumstances,  and  also  that,  the  decisions  concerning  them 
have  a  strong  Influence  In  the  discussions  about  any  other 
factors . 

In  the  determination  of  a  prior  distribution  for  the  targets 
position.  It  Is  essential  to  keep  In  mind  that  although  sub¬ 
jective  judgement  play  a  central  role  In  this  matter,  this  dees 
not  mean  that  It  Is  just  a  guessing  process.  On  the  contrary. 
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by  observing  some  rules  In  analyzing  all  the  relevant  Informa¬ 
tion  about  the  target.  It  Is  possible  to  obtain  a  distribution 
which  from  the  viewpoint  of  the  search  model,  Is  accurate 
enough. 

In  general.  It  Is  possible  to  think  of  two  basic  situa¬ 
tions  at  the  starting  point  of  a  search  operation.  Either  no 
relevant  Information  about  the  target's  last  position  Is 
available,  or  a  datum  Is  somehow  known.  In  the  first  case, 
the  uniform  distribution  Is  the  appropriate  one.  The  problem 
Is  then  to  determine  the  position  and  the  dimensions  of  the 
search  area.  Usually,  In  this  case  a  rectangle  Is  constructed 
and  adopted  as  S.A.,  for  the  sake  of  simplicity.  The  location 
and  dimensions  of  this  rectangle  will  rely  strongly  on  the 
analysis  of  the  available  Information,  and  mainly  on  expert's 
opinions.  Before  seeing  an  example.  It  Is  necessary  to  recall 
that  the  mentioned  rectangle  must  contain  the  target  with 
probability  1,  and  therefore,  the  less  Information  available, 
the  larger  the  S.A.  will  be,  and  consequently  more  effort  must 
be  placed  there.  In  order  to  assure  a  determinated  level  of 
probability  of  detection. 

EXAMPLE  1 

Consider  Figure  I,  where  "A"  and  "B"  are  airports,  and 
suppose  that  an  airplane  left  "A"  and  did  not  arrive  at  "B". 

Its  last  position  is  unknown  (from  a  very  rigorous  point  of 
view,  "A"  Is  the  aircraft  last  position,  however,  for  practical 
purposes,  this  kind  of  reasoning  Is  useless). 
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Figure  11-16 


Flight  experts  are  asked  to  Inform  what  Is  the  maximum 
error  of  navigation  for  the  aircraft  in  question.  Meteorology 
people  give  the  direction  and  Intensity  of  winds  In  the  area 
during  the  time  of  flight.  From  these  considerations  It  Is 
concluded  that  AB'  and  AB"  are  the  extreme  paths  for  the 
airplane.  Further,  the  dashed  regions  centered  at  "A"  and 
MB"  are  the  approach  control  zones  of  the  airport,  and  the 
plane  has  been  seen  leaving  the  "A"  zone,  and  did  not  enter 
the  "B"  zone. 
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Flight  control  people  assure  that  It  Is  Impossible  to 
cross  enter  zone  without  being  detected.  Finally,  considering 
that  between  MAM  and  "B",  there  is  no  place  for  an  emergency 
landing,  the  conclusion  Is  that  the  plane  crashed  somewhere, 
and  again,  by  consulting  flight  experts.  It  Is  determined  the 
maximum  distance  of  landing  from  each  point  of  the  extreme 
courses,  under  the  prevailing  conditions  of  wind.  The  search 
area  of  Figure  16  Is  then  a  result  of  all  the  above  processes. 
It  should  be  noticed  that  at  this  point,  all  the  work  of  the 
search  planner  Is  to  obtain  from  the  experts  the  necessary 
amount  of  Information  In  order  to  draw  the  minimum  search  area 
which  contains  the  target  with  probability  one. 

When  a  last  position  Is  known,  the  next  question  concerns 
how  It  was  determined.  The  reasons  Is  that  there  Is  always 
an  error  In  navigation  and  the  probability  distribution  for 
the  target's  location  will  reflect  this  error.  In  fact,  when 
a  last  position  Is  known,  the  target's  position  distribution 
can  be  viewed  as  the  error  In  navigation  distribution.  There¬ 
fore,  It  Is  by  the  knowledge  of  how  the  mentioned  last  position 
was  obtained  that  It  is  possible  to  obtain  the  form  and  the 
parameters  of  a  corresponding  probability  distribution.  In 
fact,  for  a  large  class  of  situations,  a  huge  mass  of  data 
from  past  experience  allows  a  search  planner  to  obtain  very 
accurate  results  in  fitting  a  distribution  to  a  specified 
situation.  Reference  [3]  contains  tables  and  curves  that  can 
be  used  In  many  cases,  and  which  are  results  taken  from  past 
experience.  Perhaps  these  distribution  fitting  procedures  can 
be  made  more  clear  with  the  next  examples. 
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Suppose  a  small  ooat  asks  for  help,  and  gives  her  actual 
position.  Evidently  there  Is  an  error  In  this  position,  and 
this  error  Is  assumed  normal.  Reference  [3]  tabulates.  In  Its 
chapter  6,  a  "total  error  In  position,"  denoted  "c",  which 
stands  for  the  radius  of  the  circle  which,  with  probability 
0.5,  contains  the  target.  Assuming  that  the  position  was  based 
in  dead  reckoning,  four  hours  after  last  fix  and  that  the 
average  speed  for  that  type  of  boat  Is  10  knots,  the  value 
for  the  "total  error  In  position"  is: 

C  e  46  miles 

From  the  properties  of  the  bivariate  normal  distribution, 
already  discussed  In  II. 

°x,  ■  °X2  ■  TTTJf 
or,  In  this  case: 


46 


TITS' 


38.9  miles 


Therefore,  under  the  assumption  of  normality,  If  the 
search  area  Is  centered  In  the  last  position,  the  target's 
location  distribution  Is  a  circular  normal,  with  mean  zero  and 
standard  deviation  38.9  miles.  A  search  grid  can  be  constructed 
If  necessary,  by  the  use  of  the  results  presented  In  II 

It  Is  Interesting  to  mention,  that  the  drift  effect  can 
be  computed  In  cases  like  this,  by  making 
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Where  X  Is  the  former  value  for  c  (46  miles  In  this  example), 
and  d  Is  a  "drift  error"  which  can  be  found,  for  many  types 
of  targets  and  meteorological  and  oceanographic  conditions, 

In  chapter  6  of  reference  [3]. 

Observe,  however,  that  a  drifting  target  Is  In  fact  a 
moving  target,  and  that  the  described  procedure  Is  a  practical 
approximation  that  allows  the  search  planner  to  deal  with  a 
moving  target  In  the  same  way  as  with  a  stationary  one,  just 
by  Increasing  the  variance  of  Its  location  distribution.  This 
approximation  .Is  only  valid  within  the  limits  of  the  small 
velocity  that  characterizes  drifting. 

Finally,  the  assumption  of  normality  comes  from  the 
generalized  notion  that  errors  In  navigation  are  In  fact  sums 
of  many  errors  (Instrument  reading  errors,  plotting  errors, 
etc.),  and  thus,  due  to  the  central  limit  theorem  (see  reference 
[4],  It  Is  reasonable  to  talk  of  normality. 

EXAMPLE  3 

THE  SEARCH  FOR  THE  SSN  THREASHER 
This  Is  a  real  life  example,  which  Is  discussed  in  full 
detail  In  reference  [7].  Also,  the  Interested  reader  can  find 
Important  material  about  under  water  search  In  reference  [8]. 

The  Thresher  was  lost  In  approximately  8,400  feet  of  water 
during  her  sea  trials,  after  an  overhaul,  near  the  gulf  of 
Maine.  In  establishing  her  location  distribution,  the  analysts 
In  charge,  considered  the  following  points: 


i 
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(1)  Last  Known  Position 

At  0917  on  April  10,  1967,  the  USS  Skylark  received 
a  UQC  (underwater  telephone)  message  from  Thresher  reporting 
that  she  was  In  trouble  for  maintaining  her  trim.  A  few  seconds 
after.  Skylark  heard  breakup  noises  over  the  UQC.  At  this 
same  time,  Skylark  determined  her  own  position  by  using  a  loran 
alfa.  Then,  at  this  point,  the  following  are  the  relevant 
facts . 

_  The  Skylark  position,  determined  by  the  loran  alfa, 

which  has  a  ".total  error  In  position"  (see  last  example),  of 
1.5  miles . 

_  The  nominal  range  of  UQC,  which  Is  5,000  yards. 

Therefore,  by  the  same  reasons  discussed  In  example  2,  the 
Skylarks'  position  can  be  described  by  a  circular  normal  dis¬ 
tribution  centered  at  the  loran'  fix,  and  with  standard 
deviation  given  by: 

_  _  —  1.5 

°X1  ax2  \TTT5’ 

On  the  other  hand,  by  the  knowledge  of  the  UQC  range,  It 
was  assumed  that  the  Thresher's  position  had  a  uniform  distri¬ 
bution  Inside  of  a  circle  whose  radius  Is  5.000  yards,  and 
centered  In  each  possible  Skylark's  position.  This,  leads  to 
a  compound  distribution,  as  discussed  in  II-F  for,  In  order  to 
compute  the  probability  In  any  cell  In  the  search  space.  It  Is 
necessary  to  deal  with  a  uniform  density  conditioned  by  a 
normal  one. 
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(2)  Horizontal  and  Vertical  Motion  Before  and  After 
the  Last  Known  Position 

Although  the  various  courses  steered  by  the  Thresher 
before  the  accident  was  known,  no  factual  Information  con¬ 
cerning  speeds  before,  and  courses  and  speeds  after  the 
accident  were  available.  Also  nothing  was  known  for  sure 
about  vertical  motion.  However,  by  using  experts'  opinions, 
a  maximum  range  from  the  last  position,  that  the  submarine 
could  be  carried  to  was  estimated.  Notice  that  this  maximum 
range  also  condition  the  probability  distribution. 

(3)  011  Slick 

A  diesel  oil  slick  was  noticed  nearby  the  Skylark's 
position.  However ,  no rel 1  able  Information  could  be  extracted 
from  the  experts,  because  of  the  great  uncertainty  about 
winds,  currents,  etc.,  In  the  area. 

By  using  the  Information  contained  In  (1)  and  (2),  and 
disregarding  (3),  the  search  grid  was  constructed  as  follows: 

I.  A  search  area  of  10  X  10  miles  was  established.  This 
was  done  by  computing  the  size  of  the  0.9  probability  square, 
for  the  compound  distribution  (see  II-F  example,  for  the  case 
of  a  circular  normal},  and  adding  the  range  predicted  In  (2), 
by  the  experts. 

II.  The  search  grid  was  laid  out  by  dividing  the  search 
area  In  squares,  each  one  mile  on  a  side,  and  computing  the 
respective  probabilities.  Finally,  the  following  observations 
concerning  to  the  construction  of  search  grids  are  of  Interest: 
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(1)  Whenever  possible  the  coordinate  system  for 
the  grid  should  be  the  north/south  system  of  the  mercator 
chart  for  the  are  In  question,  and  the  units  of  measure  the 
same,  l.e.,  for  the  axis  longitude,  for  the  x2  axis  latitude. 

(2)  The  choice  of  the  dimension  for  each  cell  Is 
dictated  by  the  ratio  between  the  standard  deviation  of  the 
navigation  error  (on)  and  the  sweep  width  (w).  This  Is  due  to 
the  fact  that  If  oN/W  Is  large,  then  It  Is  very  likely  that 
appreciable  amounts  of  effort  Intended  for  one  cell  will  be 
placed  outside  of  that  cell.  Reference  [7]  has  a  nomograph 
that  can  be  used  for  determining  the  correct  dimensions. 

Recall  that  values  of  on,  for  several  cases,  can  be  obtained 
by  dividing  the  "C"  factors  tabulated  In  references  [3]  and 
[7]  by  1.18.  In  reference  [3],  "C"  Is  called  "total  error  In 
position,"  and  In  reference  [7],  "circular  probable  error." 

(3)  Although  the  calculations  of  each  cell's 
probability  presented  In  II-F  might  seem,  at  first  sight,  very 
straightforward,  In  fact  they  are  usually  not.  For  example, 

In  the  search  for  the  Scorpion  (see  reference  [8],  nine 
scenarios  were  considered.  Besides  the  obvious  volume  of 
calculations.  It  Is  not  common  that  the  distribution  of  a 
random  variable,  conditioned  by  many  others,  can  be  written 

In  a  closed  analytical  form.  In  this  case,  the  correct  approach 
Is  the  use  of  a  Monte  Carlo  method.  A  computer  simulation  for 
generating  a  search  grid  Is  described  In  reference  [7].  A 
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complete  list  of  the  program,  which  Is  written  In  FORTRAN  IV, 
and  that  can  be  used  In  a  wide  range  of  situations,  can  be 
found  In  Appendix  0,  of  that  reference. 

Item  B 

After  the  target's  location  probability  distribution  Is 
known,  the  next  step  for  the  search  planner  Is  to  evaluate  what 
can  be  done  with  the  resources  he  has.  This  Is  a  very  Important 
step,  because  It  may  be  that  he  does  not  have  means  enough  to 
accomplish  his  mission.  Suppose,  for  Instance,  that  under  the 
situation  described  In  I-A  example,  we  want  a  search  with 
probability  0.95  of  finding  the  target.  As  can  be  recalled, 
this  Is  Impossible  with  the  resources  available,  which  did  not 
allow  more  than  0.45. 

Although  the  very  first  consideration  In  evaluating  a 
sensor  system,  should  be  with  the  LDF  to  be  employed,  the 
random  search  formula,  being  a  conservative  estimate  of  search 
efficiency,  should  be  employed. 

Determination  of  the  sweep  width  then  becomes  the  first 
step.  The  best  way  to  obtain  a  value  for  a  sensor's  sweep 
width  Is  to  integrate  Its  lateral  range  function,  as  described 
in  II-E.  However,  In  many  cases,  the  lateral  range  curve  Is 
not  available  and  the  following  are  special  guidelines  to 
deal  with  such  situations; 

(1)  For  visual  search,  use  the  numbers  of  the  table  In 
Figure  7-2  reference  [3]. 

(2)  For  electronic  devices,  the  following  procedures, 
listed  ir  order  of  preference,  are  recommended: 
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I.  When  the  minimum  detection  range  Is  known: 

W  *  (1.5)  (minimum  detection  range) 

II.  ^hen  average  detection  range  Is  known: 

W  «  (1.2)  (average  detection  range) 

III.  When  maximum  detection  range  Is  known: 

W  *  (0.8)  (maximum  detection  range) 

IV.  When  no  detection  range  Is  known: 

W  ■  (0.4)  (horizon  range) 

Notice,  however,  that  the  above  rules  of  thumb  should  be  used 
carefully,  and  their  results  always  checked  by  common  sense.  For 
Instance,  suppose  that  a  search  planner  wants  to  estimate  the 
sweep  width  for  an  airplane  radar,  operating  at  3,000  feet  of 
altitude  In  searching  for  a  small  boat.  By  rule  IV,  the  sweep 
width  Is  25.6  miles,  which,  for  a  small  boat  Is  exaggerated. 

The  visual  sweep  width  for  the  best  meteorological  con¬ 
ditions  for  an  airplane  at  3,000  feet,  searching  for  a  boat  of 
less  than  30  feet  long  Is  7.3  miles  (see  Table  7-2  In  reference 
[3]),  and  this  Is  a  good  estimate  of  the  radar  sweep  width  In 
this  case.  However,  If  Instead  of  a  small  boat,  the  plane  were 
searching  for  a  10,000  ton  ship,  then  the  25.6  miles  of  rule  IV 
would  be  a  good  guess. 
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Anyway,  reference  [3]  should  be  consulted  whenever  such 
rules  are  to  be  used. 

(3)  For  the  special  case  of  underwater  search,  reference 
[7]  states  very  interesting  procedures  for  estimating  sweep 
widths  of  sensors  such  as  underwater  cameras,  magnetometers, 
sonars  and  so  on.  It  should  be  mentioned  that  when  performing 
an  underwater  search,  the  sweep  width  Is  often  treated  as  a 
random  variable.  For  Instance,  consider  the  simple  case  where 
for  a  given  sensor  the  sweep  width  Is  w^  If  the  target  Is  just 
seated  In  the  bottom,  or  w2  If  the  target  Is  burled  In  the  mud. 
If  p  stands  for  the  probability  of  being  seated  In  the  bottom, 
1-p  Is  the  probability  of  being  buried  In  the  mud.  The  pro¬ 
bability  of  detection  Is: 

P(Z)  -  p (1 -EXP[1 w, Z3 )  +  (1-p)  (1-EXP(-w2Z)) 

observe  that  If  the  mean  value 

w  *  pw-j  +  (1-p)  w2 

1r  used  to  compute  P(Z),  the  expression 
P(Z)  *  1 -EXP(-wZ) 

Leads  to  a  different  result  than  (*),  which,  as  can  be  proved, 
Is  the  correct  one.  Usually,  in  underwater  search,  random 
sweep  widths  are  approximated  by  gamma  distributions.  A 
complete  discussion  can  be  found  In  reference  [2j. 
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A  point  of  Interest  Is  the  evaluation  of  the  overall 

sweep  width  when  more  than  one  sensor  Is  looking  for  the  same 

target.  In  order  to  derive  such  a  result,  ler  w  be  the  sweep 
th  J 

width  of  the  j  sensor,  and  N  Its  maximum  range.  Clearly, 

the  ratio  8j  *  }j^ 

can  be  thought  as  a  measure  of  the  probability  of  detection 
Inside  that  maximum  range.  If  n  sensors  are  used  in  a  group, 
then  the  combined  sweep  width  is  given  by: 

w  -  w,  +  w2  (i-3-j )  +  ---  wn(i-b1)  (i-e2)  ---  0 “3n_! ) 

where  W>W  >--->W 
3  ~  “  N 

A  proof  for  this  statement  can  be  found  In  reference  [7]. 

Another  relevant  point  In  determining  sweep  width  has  to  do 
with  the  targets'  physical  dimensions.  This  Is  an  entry  for 
many  tables  (table  7*2,  reference  [3],  for  example,  divides 
the  targets  In  seven  categories,  ranging  from  life  rafts  to 
ships  over  10,000  tons  of  displacement).  In  the  Thresher's 
case,  the  sweep  width  was  det  ermlned  assuming  a  268  by  300 
foot  target  which  were  the  submarine's  dimensions.  However, 
after  the  search.  It  was  concluded  that  due  to  the  Impact  on 
the  bottom  and  the  consequent  breakup,  the  light  debris 
(twisted  metal,  etc.)  scattered  by  the  submarine,  together 
with  main  parts  of  her  hull,  were  In  fact  a  target  of  1,000 
X  5,000  yards.  Therefore,  the  sweep  width  adopted  was  smaller 
than  the  real  one.  Reference  [7],  In  its  section  4.8,  discusses 
In  detail  this  point  and  tabulates  some  of  the  results.  Finally, 
an  approximated  expected  time  for  the  search  must  be  computed. 
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The  word  "approximate"  was  used  because,  at  this  point,  the 
search  planner  does  not  know  yet  what  search  plan  he  will 
employ.  In  fact,  the  choice  of  a  search  plan  will  depend  on 
this  first  rough  estimate  of  search  time.  Anyway,  by  using 
the  methods  and  formulas  presented  In  II-I,  It  Is  possible  to 
obtain  a  reasonable  approximation  to  the  expected  search  time. 
Notice  that  the  possibility  of  false  targets  must  be  carefully 
considered  at  this  point;  If  there  Is  a  considerable  density 
of  false  targets  In  the  area,  the  broad  search  time  can  be  a 
very  misleading  measure  of  time  to  be  spent. 

ITEH  C 

Suppose  that  a  search  operation  costs  an  amount  6  per 
unit  of  time,  and  that  the  value  of  the  target  Is  A.  If  T* 

Is  the  expected  time  to  find  the  target,  then  e*,  the  expected 
cost  Is: 


e*  -  eT* 

Therefore,  In  order  to  make  a  "profit",  the  search  must  begin 
only  If: 


A  >  0T* 

Observe  that  T*  can  be  computed  as  was  done  In  1 1 > I , 
either  under  the  assumption  that  the  search  continues  until 
the  target  Is  found,  or  until  a  certain  level  of  probability 
Is  reached.  In  the  latter  case,  0*  will  stand  for  the  cost 
to  reach  that  amount  of  probability.  At  this  point,  an  example 
seems  helpful. 
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EXAMPLE  4 

Suppose  that  after  an  exercise,  a  torpedo  sank.  If  Its 
value  Is  10,000  dollars,  the  cost  for  searching  Its  4,000 
dollars/hour,  and  the  available  sensor  has  a  speed  of  4  knots 
and  a  sweep  width  of  0.5  miles,  should  a  search  planner  order 
the  search  to  begin  If: 

(a)  the  torpedo's  location  Is  uniformly  distributed 
Inside  a  2  X  2  mile  search  area. 

(b)  the  torpedo's  location  has  a  circular  normal  distribu¬ 
tion  with  o*1  mile?  For  case  (a),  the  expected  time  Is: 

A  4 

T*  «  —  *  -  *2  hours 
NV  2 

therefore 

6*  -  4.000  X  2  -  8.000  dollars 

and 

V  -  10.000  >  0*  *  8.000 

Then  the  correct  decision  Is  to  start  the  search.  For  case  (b) 
the  expected  time  Is: 

2*o 

T  -  -  ■  3.14  hours 

WV 

Therefore 

6*  •  4.000  X  3.14  -  12560  dollars 
and 

V  -  10.000  <  0*  -  12560 

So  the  best  decision  Is  to  abandon  the  torpedo. 


112 


In  real  life  situations,  the  big  problem  Is  to  assign 
values  for  6  and  A.  However,  the  above  guidelines  are  useful 
whenever  such  a  decision  must  be  made. 

ITEM  D 


The  reasons  for  terminating  search  can  be  viewed  In  the 
same  context  as  the  reasons  for  beginning  search. 

In  fact.  If  the  first  Increment  of  effort  Is  spent  because 
the  expected  return  exceeds  the  expected  cost,  then  the  logical 
time  for  concluding  the  search  Is  when  the  expected  return  of 
additional  Increment  of  effort  Is  smaller  than  Ite  expected 
cost.  Again,  as  was  pointed  out  In  Item(C),  the  big  problem 
In  real  life  Is  to  estimate  value  and  costs.  However,  In  the 
stopping  case.  It  Is  possible  to  argue,  aft  er  the  probability 
of  detection  has  r  eached  a  huge  value  and  the  target  still 
has  not  been  found,  that  there  Is  statistical  evidence  to 
Indicate  that  some  assumption  Is  Incorrect  and  therefore  to 
justify  termination  of  the  search,  or  at  least  to  reevaluate 
the  Initial  assumptions  before  search  continues.  Using  the 
same  notation  as  In  Item  (c),  and  for  the  case  of  the  circular 
normal  prior,  search  must  stop  by  Ime  given  by: 


-  1)2T» 


0 

For  the  case  of  the  uniform  prior,  where  the  returns  are  con¬ 
stant,  the  rule  Is  to  stop  when  the  probability  of  detection 
reaches  some  fixed  value,  say  0.95. 
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